Mon. Not. R. Astron. Soc. OQO,[7pT](201 1) Printed 10 October 201 1 (MN KTbX style file v2.2) 



The impact of CCD radiation damage on Gaia astrometry: 

I. Image location estimation in the presence of radiation damage 

T. Prod'homme^*, B. Holl^, L. Lindegren^, and A.G.A. Brown^ 

^Leiden Observatory, Leiden University, P.O. Box 9513, 2300 RA, Leiden, The Netherlands 
^Lund Observatory, Lund University, Box 43, 22100 Lund, Sweden 

o 

Accepted 201 1 October 4. Received 201 1 April 6; in original form 201 1 April 6 

o 

^ ABSTRACT 

t — The Gaia mission has been designed to perform absolute astrometric measurements with un- 

precedented accuracy; the end-of-mission parallax standard error is required to be of the order 
' of 10 micro-arcseconds for the brightest stars (V ^ 10) and 30 micro-arcseconds for a G2V 

^5 type star of magnitude 15. These requirements set a stringent constraint on the accuracy of the 

estimation of the location of the stellar image on the CCD for each observation: e.g., 0.3 milli- 
(~| arseconds (mas) or 0.005 pixels for the same V = lb G2V star. However the Gaia CCDs will 

suffer from charge transfer inefficiency (CTI) caused by radiation damage that will degrade 
' the stellar image quality and may degrade the astrometric performance of Gaia if not properly 

^ addressed. For the first time at this level of detail, the potential impact of radiation damage on 

the performance of Gaia is investigated. In this first paper we focus on the evaluation of the 
CTI impact on the image location accuracy using a large set of CTI-free and damaged syn- 
i— ' thetic Gaia observations supported by experimental test results. We show that CTI decreases 

the stellar image signal-to-noise ratio and iiTeversibly degrades the image location estimation 
precision. As a consequence the location estimation standard errors increase by up to 6% in 
the Gaia operating conditions for a radiation damage level equivalent to the end-of-mission 
accumulated dose. We confirm that in addition the CTI-induced image distortion introduces 
ly-^ a systematic bias in the image location estimation (up to 0.05 pixels or 3 mas in the Gaia 

operating conditions). Hence a CTI mitigation procedure is critical to achieve the Gaia re- 
quirements. We present a novel approach to CTI mitigation that enables, without correction 
|— ^ of the raw data, the unbiased estimation of the image location and flux from damaged obser- 

I vations. We show that its current implementation reduces the maximum measured location 

bias for the faintest magnitude to 0.005 pixels (^4x lO^"' pixels at magnitude 15) and that the 
• • Gaia image location estimation accuracy is preserved. In a second paper we will investigate 

^ how the CTI effects and CTI mitigation scheme affect the final astrometric accuracy of Gaia 

by propagating the residual errors through the astrometric solution. 

^ Key words: instrumentation: detectors - space vehicles - astrometry - methods: numerical 

- methods: analytical - methods: statistical 
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1 INTRODUCTION 

Gaia is a European Space Agency mission that aims to create the 
most complete and accurate stereoscopic map to date of the Milky 
Way, containing parallaxes, proper motions, radial velocities, and 
astrophysical parameters for one billion stars, one percent of the 
estimated stellar population in our galaxy ( [Perryman et al.|2001a| 
[Lindegren et al.|2008| ). Due to the satellite's constant spinning mo- 
tion, the determination of the astrometric parameters ultimately 
comes down to measuring very precisely the time tobs at which a 
particular star crosses a fiducial line on the focal plane ( [Lindegren] 
|& Bastian|2011| [Bastian & Biermann|2005l l. The required astro- 

* E-mail: prodhomme@strw.leidenuniv.nl 



metric precision is extreme, e.g., the end of mission parallax un- 
certainty for a star of magnitude V = 15 is required to be better 
than 25 micro-arcseconds (/ias|^ In order to determine tahs, one 
needs to measure the image location on the Charge-Coupled De- 
vice (CCD) relative to the instrument axes. As a consequence the 
required astrometric accuracy sets a direct and stringent require- 
ment on the residual image location uncertainty per CCD star tran- 
sit. In the left part of Table [T| we detail the end-of-mission paral- 
lax standard error, o.^, as function of stellar magnitude and typ^ 

^ A list of acronyms is provided in Table|2]at the end of the paper 
^ Updated estimates of the science performance are given on: 
http : / / www . rssd . esa . int/ index . php ?project= | 
GAIA&page=Science_Perf ormancej 
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computed using |de Bruijne| ( (2009[ >. These predicted standard errors 
do include the increased photon noise due to the radiation dam- 
age induced charge loss, but not the residual bias-calibration er- 
rors considered in the present paper (except for a general contin- 
gency margin of 20%). For this paper we are interested in the mean 
image location uncertainty per CCD star transit a^. that would be 
needed to reach a given targeted parallax accuracy. Based onjde] 
|Bmyne]|2005 1 we estimate the corresponding 'requirement' on the 
image location uncertainty shown in the right part of Table [T] and 
computed as: CTk = \/Nahs o:^ /{m g^), with A^obs the average 
number of astrometric observations per star (662), m the end-of- 
mission scientific contingency margin which is 1.2, and g^j the ge- 
ometrical parallax factor which is 2.08 for the Gaia solar aspect 
angle ^ — 45°. This formula has also been used to compute the 're- 
quirement' curve as function of Gaia G-bancj^magnitude shown in 
several figures throughout this paper No spectral type distinction is 
needed when these uncertainties are expressed in G because they 
are virtually independent of spectral type. Note that the computed 
location uncertainties do not contain the 20% contingency margin, 
making them very stringent. 

During the 5 year mission life time, solar wind protons will 
collide with Gaia's focal plane and create electron traps in the 
CCDs by displacement damage. These radiation-induced traps 
drastically increase the CCD charge transfer inefficiency (CTI) and 
will lead to a significant loss of signal for all Gaia measurements 
by stochastically capturing and releasing signal electrons. The re- 
sulting electron redistribution will also distort each stellar image. 
The CTI effects are expected to significantly contribute to the er- 
ror budget of all the Gaia measurements (astrometric, photometric, 
and spectroscopic), especially if not properly taken into account in 
the data processing. 

We present here the first part of a detailed evaluation of the 
impact of the radiation damage effects on the final accuracy of the 
Gaia astrometric measurements. This paper focuses on the effect 
of CTI on the image location accuracy. Studying the accuracy of 
a measurement is a rather complex enterprise. Hence we present 
in the following section the overall applied methodology and the 
different steps of this study. In a second paper, we will investigate 
how the final Gaia astrometric accuracy is affected by the CTI- 
induced errors at the image processing level as characterized in this 
study. 



2 OVERALL METHODOLOGY 

The use of synthetic data To evaluate the impact of the CTI ef- 
fects on the image location accuracy, we apply the Gaia image pa- 
rameter estimation procedure (Section|4j to a large data set of simu- 
lated CTI-free and CTI-affected observations (from here on the lat- 
ter are referred to as 'damaged observations'). The use of synthetic 
data presents several fundamental advantages compared to the use 
of experimental data: while in experimental studies, the true image 
parameters, the instrument model or Point Spread Function (PSF), 
and the different noise contributions need to be estimated, in the 
simulation, these are known parameters. Hence the uncertainties 
related to the estimation of such parameters cannot bias the result 



^ The Gaia G-band magnitude is a broad-band, white-light magnitude in 
the wavelength range 300 - 1000 nm detined by the telescope transmission 
and CCD quantum efficiency. G = V for an un-reddened AOV star ^Jordi| 
|et al.|2010[|Perryman et al.|2001b| . 



Parallax accuracy target Con'esponding CCD image 
standard error cr^ [/^as] location uncertainty [mas] 



Type 


BIV 


G2V 


M6V 


BIV 


G2V 


M6V 


V-G 


0.03 


0.16 


2.18 


0.03 


0.16 


2.18 


V = W 


6.9 


7.0 


7.3 


0.072 


0.073 


0.076 


V = 15 


25 


24 


10 


0.26 


0.25 


0.11 


V = 20 


322 


300 


102 


3.3 


3.1 


1.1 



Table 1. On the left side we tabulate the Gaia target performance of the sky 
averaged end-of-mission pai'allax standard error, (Tu, (in /las), as function of 
spectral type and Johnson V-band magnitude (see |de Bruijne|2009) . V—G 
allows for the conversion from V to the Gaia G-band magnitude. On the 
right we give the corresponding mean CCD image location uncertainty, cr^ 
(in mas), that would result in the parallax standard en'or on the left when 
observed 662 times (the average number of astrometric observations per 
star). See Section^for more details. 

of our study. Furthermore only simulation can allow the determina- 
tion of the absolute image location bias and the associated standard 
errors as this requires the knowledge of the true image location. 
Finally, by using synthetic data one can compute the statistical un- 
certainties on the measured image location bias and standard errors, 
and this at any precision level just by increasing the number of sim- 
ulated observations for a particular set of conditions. In Section|3] 
we detail the simulation of Gaia-like observations in the absence 
and presence of radiation damage for different stellar magnitudes, 
image widths, background and readout noise levels. 

The bias allows for the quantification of the trueness of an es- 
timation, and the standard errors for the quantification of the esti- 
mator precision. If an estimator delivers bias-free estimations, then 
its standard errors can also be regarded as a means to quantify the 
accuracy of this estimation. In the following we will thus make the 
important distinction between precision and accuracy when it is 
justified. 

The Gaia image parameter estimation procedure In the Gaia 
data processing, the image location and flux are estimated or 'self- 
calibrated' through the use of an iterative procedure, that allows for 
the successive determination and improvement of the PSF and the 
image parameters without prior knowledge. A detailed description 
of this procedure is provided in Section|4] In Sections [5. 3| and [54] 
the procedure is applied to the data set of simulated CTI-free ob- 
servations, in order to verify that, in the absence of CTI, the Gaia 
image parameter estimation procedure performs efficiently, accord- 
ing to expectations. 

The Cramer-Rao bound Assessing the efficiency of the Gaia im- 
age parameter estimation procedure necessitates the computation 
of the theoretical limit to the image location accuracy in the Gaia 
observing conditions. This theoretical limit corresponds to the ulti- 
mate accuracy achievable by any bias-free estimator. It is set by the 
Cramer-Rao bound, described in Section [5T| We thus compute the 
Cramer-Rao bound as a function of magnitude (G), image width, 
background, and readout noise level (Section [5^ and subsequently 
use it as a reference in a comparison with the standard errors of 
the estimated image parameters. The Cramer-Rao bound is also re- 
quired to assess the impact of the CTI effects independently from 
any estimation procedure. 

The radiation damage impact on the image location estimation 
In Section [53] we use the set of damaged observations to demon- 
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strate that the image distortion and the charge loss induced by CTI 
imply an irrevocable loss of accuracy in the image location deter- 
mination. This loss of accuracy, which directly affects the perfor- 
mance budget of Gaia, is independent from any estimation method 
and can only be avoided by physically preventing charge trapping. 
This is done by optimizing the hardware (e.g., the CCD operating 
temperature) and using hardware countermeasures such as the pe- 
riodic injection of artificial charges, or the use of a supplementary 
buried channel, an extra doping implant in each pixel that confines 
small charge packets. Taking into account these countermeasures in 
the simulation of the damaged observations allows us to verify that 
they indeed substantially contribute to diminish the CTI effects on 
the theoretical image location accuracy. However, Sections [5.6| and 
|5.8| show that it is not possible to rely solely on these countermea- 
sures: indeed we find that if the CTI effects are not properly taken 
into account in the image parameter determination the image loca- 
tion bias can be as large as 10 mas for a star of magnitude 15 (to 
be compared to the requirement of ~0.08 mas, see Table[TJ. This is 
in agreement with the experimental tests performed on Gaia irradi- 
ated CCDs (see Section [577} and confirms that the CTI effects need 
to be addressed by the Gaia data processing in order to achieve the 
mission requirements. 

Mitigating the CTI effects The software mitigation of CTI effects 
is a complicated task. Several schemes have been discussed in the 
literature; they usually imply the direct correction of the raw data in 
the context of photometry-based measurements. In Section [6T[ we 
review the different potential schemes, and where they intervene 
in the data processing chain. Then, in Section [6!2l we present and 
motivate a novel approach to CTI mitigation that does not involve 
a direct correction, so that the noise properties of the raw obser- 
vations remain unchanged. This approach, developed by the Gaia 
Data Processing and Analysis Consortium (DPAC, |Mignard et al.| 
|2008[ l, relies on the forward modelling of each observation includ- 
ing the CTI distortion, such that the true image parameters can be 
directly estimated from the damaged observations. The modelling 
of the distortion of the stellar image is performed thanks to a fast 
analytical CTI effects model, a so-called Charge Distortion Model 
(CDM). The success of this CTI mitigation approach depends on 
the performance of such a model. In Section [63l the potential ac- 
curacy of this approach is assessed in the case of an ideally cali- 
brated CDM. Finally, using the current best CDM candidate ( |Short| 
|et al.|[20T0{ [Prod'homme et al.||2010[ and Section [6!4](, we apply 
the image parameter estimation procedure and the DPAC CTI mi- 
tigation scheme to the set of damaged observations and show that 
one can recover the CTI-induced image location and flux bias (see 
Sections |6.6[ Only then are we able to answer the question: 
does the current Gaia image location procedure combined with the 
presented CTI mitigation scheme allow an unbiased estimation of 
the image location with a sufficient accuracy in the presence of ra- 
diation damage? 

Methodology summary 

(i) Generation of CTI-free and damaged Ga;a-like observations. 

(ii) Determination of the theoretical limit to the image location accu- 
racy in the absence of CTI by computing the Cramer-Rao bound. 

(iii) Performance assessment of the Gaia image parameter estimation 
procedure in the absence of CTI. 

(iv) Evaluation of the intrinsic loss of accuracy in the image location 



estimation induced by radiation damage by computing the Cramer- 
Rao bound for a damaged LSF. 

(v) Characterization of the CTI effects on the image parameter esti- 
mation procedure. 

(vi) Performance assessment of the Gaia image parameter estimation 
procedure in the presence of CTI and including a forward model- 
ling approach to mitigate the CTI effects. 



3 GENERATING GAM -LIKE OBSERVATIONS 

In the following, we first describe the main principles of the Gaia 
observations, and explain how we generate GaM-like reference im- 
ages. These images are used to simulate thousands of observations 
for different stellar magnitudes and operating conditions, and thus 
constitute the basis of our study. Then we detail how we simulate 
the stellar transits over a CCD. To achieve a high level of realism, 
we use a physically motivated Monte Carlo model that simulates 
the CCD charge collection and transfer, as well as the trapping pro- 
cesses, at the pixel-electrode level ( [Prod'homme et al.|2011^ . Fi- 
nally we summarize the expected CTI effects on the stellar images 
in the Gaia operating conditions and explain the choices we made 
regarding the radiation damage parameters of the simulation (trap 
species, level of radiation). 

3.1 How Gaia observes 

The Gaia spacecraft will orbit around the second Lagrange point 
(L2) and constantly spin around its own axis such that its two tele- 
scopes scan a great circle on the sky several times a day. The pre- 
cession of the spin axis changes the orientation of the consecutive 
great circles, allowing for the coverage of the whole sky in about 
six months. The measurements are recorded in a single focal plane 
consisting of 106 CCDs. Due to the satellite spinning motion, the 
star projections will not remain stationary during an observation 
but will transit the focal plane in the along-scan (AL) direction. 
The orthogonal direction is called across-scan (AC). To integrate 
the stellar images along the star transits, the CCDs are operated 
in time-delayed integration (TDI) mode. In this mode the CCD is 
constantly readout and the satellite scanning rate (and induced light 
source motion) has been synchronized with the charge transfer pe- 
riod, so that the charge profile continues to build up as the image 
travels across the CCD avoiding as much as possible image smear- 
ing. The charge transfer period is 0.9892 ms and the integration 
time 4.4 s. The observing principle of Gaia relies on differential 
positional measurements among the stars simultaneously visible in 
the two superposed fields of view. In particular the differential mea- 
surements between the two fields of view (covering arcs of about 
106.5° on the sky) are essential for the construction of a global 
reference frame and the determination of absolute parallaxes. For 
these measurements the AC component of the differential positions 
is largely degenerate with the instrument pointing, and mainly the 
AL component matters l |Lindegren & Bastian|20I ij l. Gaia is there- 
fore primarily optimized for AL measurements and the image lo- 
cation accuracy in that direction is the most critical one for the per- 
formance. 

Because of limitations on the amount of data that can be sent 
to ground only a small window around each source is read out. 
For sources fainter than magnitude 13 the windowing scheme is 
simple. In the AL direction, the window size is 12 pixels for stars 
brighter than magnitude 16, and 6 pixels for fainter stars. In the AC 
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Figure 1. One-dimensional Gaia-hke reference images generated from a 
symmetrized subset of realistic Gflfci-like line spread functions (LSFs). 
These reference images are used to reflect the range of narrow, typical, and 
wide profiles that results from different stellar types and wave front errors. 
The bottom figure shows the same profiles as the top figure for a wider AL 
coordinate range and with a logarithmic ordinate scale. 



direction the size of the on board readout windows is 12 pixels (the 
pixel size is 58.9 mas AL and 176.8 mas AC). The observations are 
then binned across-scan before being sent to the ground. Hence, for 
a particular star, a Gaia observation results in a one-dimensional 
along-scan set of electron counts that corresponds to the sampling 
of a one-dimensional point spread function or Line Spread Function 
(LSF). Note that Gaia will observe sources as bright as G = 5.7, 
but due to the relatively low number of stars between magnitude 
5.7 and 13 (~1% of the expected lO'' sources) and because of the 
use of a complicated gating scheme (to avoid pixel saturation) these 
magnitudes are ignored in this study. 

3.2 Construcing a Gaia-like reference image 

In this study we only analyze observations of point sources be- 
cause these will be used in the Astrometric Global Iterative Solu- 
tion (AGIS, [Lindegren et al.|20rT| > to calibrate the instrument and 
satellite attitude, which are then used to estimate the astrometric 
parameters for all sources. The PSF, the actual two-dimensional 
flux distribution of an unresolved star that illuminates a CCD, de- 
pends on the spectral energy distribution of the star, the CCD prop- 
erties, and the optics and its associated wavefront errors. Because 
the spectral energy distribution depends on the type of star, and the 
wavefront errors depend on the position in the focal plane, we can- 
not construct a single PSF that would be representative of the whole 
range of possible profile shapes. 

To construct a set of PSFs that samples the range of possible 
profile shapes we make use of a study by |Lindegren| l |2009) l in which 



a set of 20,000 one-dimensional line spread functions (LSFs) were 
generated that are representative of the Gaia optics and wave front 
errors, stellar spectral energy distributions, and CCD effects (e.g., 
smearing due to TDI). A set of basis functions was extracted from 
this dataset using a Principle Component Analysis (PCA) in order 
to describe the data with a minimum number of parameters. A set 
of ten basis functions was found to be sufficient to represent any 
of the LSFs with an RMS error of 10^*. To be readily usable the 
basis functions were subsequently fitted by a special quartic splin^ 
that is flexible enough to fit the data while being smooth enough to 
avoid overfitting that could result in sub-pixel position bias when 
the function is used for location estimation on pixel sampled data. 

For our study we decided to use three reference images. First 
of all we construct a reference profile based on a selection of all 
profiles that have a FWHM within ± 1% of the mean FWHM 
of 1.958 pixels. These profiles are subsequently symmetrized (i.e. 
Laymm = {L{x) + L{—x))/2), avcragcd, and then fitted with the 
first four even basis functions, resulting in a symmetric mean profile 
with FWHM 1.957 pixels, from hereon referred to as the 'typical' 
reference profile. Four components are used to get a profile that is 
sufficiently close to the target mean FWHM. To represent the ex- 
tremes we introduce a 'narrow' and a 'wide' reference profile. The 
narrow and wide profiles are constructed in exactly the same way as 
the typical reference profile, only differing in the selection of LSF 
samples, which are: 90 ± 1% and 110 ± 1% of the mean FWHM 
respectively. This results in a FWHM of 1 .767 pixels for the narrow, 
and 2.161 pixels for the wide reference profile. All three reference 
images are normalized as shown in Fig.[T] 

Although the two-dimensional PSF in the Gaia focal plane 
is wider in the across-scan direction than in the along-scan di- 
rection, the pixels are shaped such that in pixel units the PSF is 
nearly identical in both directions. Because we are only interested 
in how the pixels are illuminated we can therefore construct the 
two-dimensional reference image from simply multiplying the one- 
dimensional reference image L{x) in two dimensions: 



P{x, y) ^ L{x - k) X L{y - fi) 



(1) 



In all our further analyses we assume fi — 0. Because we defined 
the zero-point of the symmetric profile L to be at the symmetry 
point, this ineans that the PSF is always in the centre of the window 
in the across-scan direction. 



3.3 Monte-Carlo simulations of observations 

The approach we have chosen for this study is to simulate a fully 
synthetic dataset using a detailed physical simulation of the photo- 
electron collection and transfer in CCDs at the pixel-electrode level 
( [Prod'honmie et al.||2011t, availabl e through the CEMGA soft- 

The model also allows for a 
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detailed treatment of radiation induced traps that capture and re- 
lease electrons and thereby distort the charge profile transferred 

^ The special quartic spline used to represent the LSF can be described 
as the convolution of an ordinary cubic spline, defined on a regular knot 
sequence with a knot separation of half a pixel, with a rectangular function 
of width equal to one pixel. This spline has the property that the sum of 
points sampled at one pixel separation is independent of the sub-pixel phase 
of the sampled points. Any 'effective' LSF, being the result of an optical 
LSF convolved with the pixel response function jAnderson & King|200Q^ , 
should have this property. Details of the special quartic spline are found in 
a technical note by .Lindegren, (,2003j. 

^ |Www . st rw . le idenuniv . nl/ "prodhomme/cemga . php| 
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through the CCD (see Section[34}. The observations are simulated 
in two dimensions: 4494 pixels in along-scan and 12 pixels in the 
across-scan direction. In the software we illuminate the CCD with a 



two-dimensional reference image P {x,y) described in Section 3.2 



The normalized reference image is scaled to produce an illumina- 
tion that corresponds to a particular stellar magnitude. The photon 
detection is modelled as a Poisson process: at each transfer step, 
the photo-electrons are generated using a random generator with a 
Poisson distribution and a mean equal to the expected number of 
collected photons (given by the reference image) within the inte- 
gration time (~lms x 4494 pixels in the case of Gain). Note that 
we control the exact along-scan location of the reference image, 
therefore allowing us to determine the exact error when a location 
estimate from the observation has been made. In the simulations we 
can optionally include a constant background. The electron packet 
transfer in the readout register is not simulated. 

The raw two-dimensional observation counts are cropped in 
along-scan direction to 80 pixels (centred around the signal) and the 
resulting 80 x 12 pixels are stored. All used reference images are 
zero for > 20, therefore any relevant signal is always contained 
in the cropped raw observation data. 

When processing an observation we load the raw two- 
dimensional pixel counts. Depending on the windowing scheme for 
this particular magnitude and CCD we crop the data around the sig- 
nal to the relevant window size and optionally bin the pixel counts 
in the across-scan direction resulting in a one-dimensional sample 
of the transit photo-electron counts {A'^fe}. Readout noise can be 
added to the counts using a normally distributed random genera- 
tor Normal(0, r^): having zero mean and standard deviation r (the 
readout noise value). 

Our total synthetic observational dataset consists of: 



(i) 3 different CCD states: CTI-free, damaged with 1 trap pixeP 
and damaged with 4 traps pixel~^ (see Section 3.4 for details about 
the damaged cases), 

(ii) 3 different reference images: narrow, typical and wide (see Sec- 
tion|T2ll, 



(iii) 2 different levels of sky background: and 0.44698 e pixel s 
(the latter corresponding to the average sky surface brightness), 

(iv) 9 different magnitudes: G= 13.3, 14.15, 15.0, 15.875, 16.75, 
17.625, 18.5, 19.25, 20.0, 

(v) the (two-dimensional) photo-electron counts of 250 CCD transits, 
each with the reference image incrementially shifted by 1/250* of 
a pixel in the along-scan direction. 

In almost all of the processing we select a unique combination of 
(i), (ii), (iii) and magnitude (iv), containing all 250 transits. The 
selection of all across-scan binned transits for a given magnitude is 
denoted as: {{Nk}}G- 

3.4 Simulation of tlie CTI effects 

At L2 the radiation environment is dominated by energetic protons 
emitted during solar flares. The proton fluence is thus governed by 
the cyclic activity of the Sun which is usually monitored through 
sunspot counts. According to the latest predictions (see e.g, |SID(>| 
|team||20m , the next peak of activity will occur during 2013 co- 
inciding with the launch of Gaia. Using the JPL 1991 model, the 
reference interplanetary proton fluence model by |Feynman et al^ 
( |I993) , taking into account the satellite design, and assuming 4 
years of operation during the solar maximum (and one year dur- 
ing minimum), the average accumulated radiation dose received by 



a CCD of the astrometric instrument is predicted to be ~ 3x 10^ 
(10 MeV equivalent) protons cm~'^ . These protons will collide with 
and displace atoms in the Gaia CCD silicon lattice, and lead to 
the creation of electron traps. These traps stochastically capture 
and release the electrons transferred in the CCD. For more infor- 
mation concerning the trapping processes see |Prod'homme et al.| 
l |2011j l and references therein. The traps originate from different 
chemical complexes generally referred to as a trap species: a sum- 
mary of the expected trap species in the Gaia CCDs is provided 
by |Seabroke et al.| ( [2008| l and |Hopkinson et aLl ( |2005| . One usually 
distinguishes between trap species with short and long release time 
constants relative to the characteristic trap-electron interaction time 
(~1 ms for the Gaia CCDs), as they have different effects on the 
measurements. The traps with short release time constants capture 
electrons from the image leading edge and redistribute them within 
the telemetry window, which induces a distortion of the charge pro- 
file. The traps with longer release time constants capture electrons 
from the stellar profile and release them outside the telemetry win- 
dow, which implies a charge loss that reduces the signal to noise 
ratio (see Fig.|2]l. The Gaia CCDs comprise two hardware CTI mi- 
tigation tools: a charge injection (CI) structure and a supplementary 
buried channel (SBC). The CI structure is located all along the first 
CCD pixel row; it is composed by a diode capable of generating 
artificial charges and a gate that controls the number of electrons 
to be injected in the first pixel row and subsequently transferred 
across the whole CCD. Charge injections temporarily fill a large 
fraction of the traps present in the CCD and effectively prevent the 
trapping of the following generated and transfered photo-electrons. 
The SBC is a second and narrower doping implant on top of the 
buried channel. It creates a deeper potential that disappears into the 
shallower but wider buried channel for charge packets larger than 
1500 e~. By concentrating the electron distribution into a smaller 
volume it minimizes the electron-trap interactions in the rest of the 
pixel volume, effectively reducing the fraction of trapped electrons 
at low signal levels (^ 1500 e^or G > 15). 

In order to obtain representative results from our study, it is 
critical to achieve a high level of realism in the simulation of the 
CTI effects on each observations. This is why we make use of 
the most detailed CTI effects model to date (Prod' homme et aT] 
|20I I[ l verified against experimental tests performed on Gaia irradi- 
ated CCDs. At each transfer step, this model simulates the capture 
and release of electrons by computing for each trap the capture and 
release probabilities according to the trap characteristics and the lo- 
cal electron density distribution taking into account the Gaia pixel 
architecture and in particular the presence of the SBC. In the fol- 
lowing we detail the considerations that led us to choose to simulate 
a unique trap species and two different radiation levels. 

During the mission CI will be performed at periodic inter- 
vals. This means that most of the traps with release time constants 
greater than the injection period will be permanently filled, as only 
a very small fraction of them will have the time to release an elec- 
tron. The current most likely value to be selected for the injection 
period is 1 s. If one neglects the serial CTI effects (occurring dur- 
ing the charge transfer in the CCD readout register), and take for 
reference the trap species as presented in |Seabroke et al.| ( |2008[ ), 
the only trap species that remains significantly active corresponds 
to the so-called 'unknown' with a release time constant r ~ 90 
ms at the Gaia operational temperature and a capture cross-section 
a = 5 X 10^^" m^. We thus decided to generate the damaged ob- 
servations using a virtual irradiated CCD containing a unique trap 
species, with these parameters. Note that the release and capture 
time constants vary exponentially with the temperature. The tem- 
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Figure 2. The two plots show the simulated damaged observation resulting from multiple transits over an irradiated CCD containing 4 traps pixel"^ at different 
magnitudes with no background (left) and a background set to the sky brightness (right). In both cases the typical reference image is used and no readout 
noise has been added. The transfer direction is from left to right. The ordinate scale corresponds to the photo-electron counts normalized by the total expected 
CTI-free flux at each magnitude. Only 25 observations per magnitude (out of 250) are shown. The reference image (line) enables us to appreciate the variation 
of amplitude in the CTI induced distortion and charge loss as a function of the signal level. Also the overall profile centroid shift is clearly visible. Note that, 
as one can observe from the left plot, a radiation dose leading to 4 traps pixel" ^ translates into a severe distortion of the Gaia measurements, pailicularly 
for fainter stars. The expected end-of-mission dose, based on the latest prediction of the next solar cycle, is however three quarters of the one we chose to 
simulate. Moreover the right plot shows that even a low level of background (~ 1 e ^ pixel^ ^ ) strongly mitigates the CTI effects by filling the traps prior to the 
star transit. 



perature over the entire Gaia focal plane is expected to deviate at 
most 5 K from the nominal operating temperature. This means that 
for different CCDs the effect of a single trap species will be dif- 
ferent. However the temperature variation over a single device is 
expected to be negligible, hence our assumption regarding a single 
trap species with a unique release time constant still holds. 

The traps filled by the CI, release their electrons and induce a 
characteristic 'release' trail after the CI. This trail changes the uni- 
form nature of the background and must be carefully taken into ac- 
count in the background estimation procedure. In order to prevent 
the CI background estimation and removal from affecting the re- 
sults of our analysis, no CI was performed before the stellar transit 
during the simulations. Hence, the trap density has to be carefully 
selected to reproduce the trap density as it is perceived by a star 
after a CI delay (or time since last CI) comparable to the CI period. 
In this way, the simulated amplitude of the CTI effects corresponds 
to the one observed in the experimental tests with CI performed in 
similar conditions. 

In a series of four different campaigns of experimental tests 
carried out on irradiated Gaia CCDs, the prime contractor for Gaia, 
EADS Astrium, investigated the performance of potential hardware 
mitigation tools and characterized the trend and amplitude of CTI 
effects on Gaw-like measurements. These campaigns are referred 
to as radiation campaigns (RC). The RCs were performed in simu- 
lated Gaia operating conditions: a CCD operated in TDI mode at a 
temperature of 163 K with a low level of background light. The de- 
vices were irradiated at room temperature with a radiation dose of 
4 X 10^ protons cm~^ (10 MeV equivalent) that corresponds to an 
upper limit to the predicted Gaia end-of-life accumulated radiation 
dose. A trap density of 4 traps pixeP^ is necessary to reproduce 
the amplitude of the CTI effects, in particular the fractional charge 
loss as observed in the second RC (RC2) from first pixel response 
measurements ( [Prod'homme et al.|2011^ . This test was performed 
with charge injections occurring every ~ 27 s. The relative image 



location bias was measured in similar conditions during the same 
campaign for a star transit occurring 1 and 27 s after the last CI. 
These results are summarized in Fig.[T2]along with the absolute lo- 
cation bias computed in this study. For a CI delay of 1 s the location 
bias is clearly smaller than for a longer CI delays (e.g., CI period 
of ~27 s), this is due to the fact that shorter CI periods maintain a 
larger portion of the traps constantly filled. As a consequence, our 
simulations were performed for two different active trap densities 
(or level of radiation damage), 1 and 4 traps pixel^^. By active we 
mean empty before the transit of the star of interest over the CCD. 
These densities reproduce the amplitude of the CTI effects as ob- 
served for short and long CI periods in the experimental tests. 

Figure [2] shows the resulting simulated CTI-induced distor- 
tion and charge loss by comparing, for different illumination lev- 
els, the simulated CTI-free observations and damaged observations 
(4 traps pixeP^) after a normalization. Note that for a unique dam- 
aged CCD containing a single trap species with fixed parameters, 
the distortion varies significantly from one signal level to the other 
and not linearly. This is due in particular to the SBC, which miti- 
gates the CTI effects at low signal level. 



4 THE GAIA IMAGE LOCATION ESTIMATION 
PROCEDURE 

4.1 Observation model: scene 

To model the flux distribution that illuminates a CCD we need 
a model of the instrument response to a point-like source, and a 
model of the actual distribution of (point) sources on the sky. The 
former has already been parameterized in Section |3.2[ it is given 
by the line spread function L when considering one dimension, or 
the PSF when considering two dimensions. Because we will mainly 
deal with one-dimensional data in this study we will hereafter only 
refer to the one-dimensional LSF For the purpose of this study 
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a simple observation model is sufficient: E{Nk) is the expected 
number of photo-electron counts in pixel k, Xk is the modelled 
photo-electron count given by a flat background /3 plus a single 
point source with flux a at location k: 

E(Nk) = Xk= aL{k~K) + l3 (2) 

Here a, k, and /? are called the scene- or image parameters. 

4.2 Maximum-likelihood estimation of ttie image parameters 

In the Gaia data processing the image location k and flux a 
will be estimated by fitting the modelled photo-electron counts 
{Xk} (Eq. [2]l to the observed photo-electron counts {Nk} using 
a Maximum-Likelihood (ML) algorithm, and this for each obser- 
vation.The image backgroimd 13 is not determined using the ML al- 
gorithm but beforehand by a more adequate method. This method 
makes use of empty telemetry windows to estimate separately the 
different main components of the background: astrophysical back- 
ground (zodiacal light, faint stars and galaxies i.e. G > 20), CI 
trails, and the CCD electronic offset. In this study we always as- 
sume that /3 is known. Another parameter that is considered to be 
known beforehand is the CCD readout noise r. Therefore, when 
estimating any of the image parameters, the true values of /3 and 
r will be used. The ML algorithm is comprehensively described in 
|Lindegren| ( [2008[ >, therefore only the main assumptions and equa- 
tions are detailed in this paper. 

According to the ML principle, the best estimate of the param- 
eter vector (here Oi — k and 62 ~ a) maximizes the likelihood 
function or equivalently the log-likelihood function: 



£{e\{Nk}) = J2^np{Nk\Xk,r) 



(3) 



with p the probability density function of the sample value given 
the modelled count and readout noise. We hence need to adopt a 
probabiMty model for the sample values. To do so we assume (i) 
that the noise is not correlated from a sample to another (already 
implicit in the sum in Eq.|3j; (ii) that the variance of the noise is 
E[{Nk — Xk)^] ~ Afc and (iii) that the sample value including 
the readout noise can be modelled as P oissonian ra ndom variables, 
Nk ~ Poisson(Afc + r^) — (see Lindegren 



2008 I. That the 



Poisson distribution is discrete, while A'^^ (obtained by correcting 
the digitized values for bias and gain) are in general non-integer, is 
not a problem as long as A'^^ + ^ 0. The continuous probability 
density function derived from the Poisson distribution is: 



p{Nk\Xk,r) = const x 



r{Nk 

and Eq.|3]can then be re-written: 

£{0\{Nk})=const+Y,[{Nk + 



-r2 + 1 



r^)ln{Xk{6)+r^) 



(4) 



-Afc(0)] (5) 



which is maximized by solving the following system of equations: 



dei0\{Nk}) 
de 



E 



jVfc-Afc(6>) dXk 
Afc(6l) + r2 de 



(6) 



These equations are non-linear and must be solved by iteration. 
Given an initial estimate 0^'^\ the linear system to be solved in 
iteration m is: 



1 (m) A o(m) 



AO'-"'' = S 



(m) 



whereupon 



j(m + l) 



j(m) 



A0 



(m) 



(7) 



(8) 



A is a symmetric positive definite matrix computed from the ex- 
pectation of the Hessian matrix; its elements are: 



- E 



1 dXk dXk 



and 



^ _ ■sp 1 dXk 



(9) 



(10) 



The iterations converge quickly if the initial estimate is reasonably 
close to the ML solution. 



4.3 First image parameter estimates and LSF model 

The ideal image model L (the true underlying flux distribution for 
each observation) corresponds to the reference image that is used to 
generate the data. During the mission, L will not be known. There- 
fore we have to estimate an image model L using the observations 
themselves. This estimation is an iterative process (Section \4.4\ , 
and successive iterations are denoted with the superscript (n) for 
n = 0, 1, . . . (not to be confused with the iterations in Eqs.|7]and[8]l. 

Given a set of transits for a certain reference image, back- 
ground and G, denoted by {{Nk}}G, how do we make the first 
estimate of the image paramete rs an d generate the first LSF model, 
L*^"' ? As mentioned in Section ' 



4.2 



the background /3 and readout 
noise r are assumed to be known already. The most straightforward 
initial flux estimate a^^^ can be made by simply taking the sum of 
the observed counts after subtracting the background. The initial es- 
timate for the image location K^°^ is determined using Tukey's Bi- 
weight centroiding algorithm (P ress et al.|1992[|Linde gren 2006). 

To generate the first estimate L'-"' we use the initial location 
estimates {k^'^^}g to relatively align the photo-electron counts of 
all selected profiles and create an oversampled profile. The creation 
of the oversampled profile is possible because each count results 
from the sampling of the reference image at a different sub-pixel 
position (Section|33J. After a background subtraction the oversam- 
pled profile is fitted by the special quartic spline to obtain Z*^"' . This 
profile estimation procedure is illustrated in Fig. [3] 



4.4 Iterative image parameter and LSF model improvement 

Once the first image model Z'"' is available, an improved estimate 
of the image parameters of each individual transit can be made us- 
ing the ML algorithm (described in Section [4~2) . This is illustrated 
in Fig.|4] Based on these improved image parameters an improved 
image model can be constructed, leading to the iterative scheme 
shown in Fig. |5] where the image parameters and image model are 
improved one after the other. Note that in the whole procedure we 
have not used any prior knowledge: everything is estimated from 
the observed photo-electron counts (i.e. 'self-calibrating'). 

After each iteration the residuals between the modelled and 
the observed photo-electron counts are monitored through the com- 
putation of the '■ 



2 

Xg 



T-lK-l 2 

— With (Ttk = Ntk + r (11) 



where T is the total number of transit profiles for a certain G, and 
K the number of along-scan pixels in each profile. For transit t 
and pixel k, Xtk and Ntk are the predicted and observed photo- 
electron counts respectively. The uncertainty atk is considered to 
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Figure 3. Diagram of the construction of the estimated LSF model L. For 
a particular reference image, background level and G, the photo-electron 
counts of all selected profiles are aligned relative to each other using the 
current best estimates of the scene parameters to create an oversampled pro- 
file, which is then fitted by the special quartic spline using a Least Squares 
algorithm, after removal of the background. 
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Figure 4. Diagram summarizing the estimation of the scene (or image) pa- 
rameters from a single observation. Modelled counts, computed using the 
latest LSF model, are fitted to the observed counts using our Maximum 
Likelihood (ML) algorithm (see Section [43) . Note that the background is 
not determined by the ML algorithm but by a dedicated procedure not de- 
tailed in this paper. In our study, the background /3 is assumed to be known. 
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iterative paratiteter update 

Figure 5. Resulting top level diagram of the Gaia image parameter esti- 
mation iterative procedure. The LSF model, L, and the scene parameter 
estimates, k, a (respectively the image location and flux), are iteratively im- 
proved by fitting the modelled counts {Aj.} to the observed counts {A^j;}. 
The modelled counts are predicted by our observation model (Eq.|2| based 
on the cuiTent best scene parameter estimates. In our study the background 
P is assumed to be known with a high accuracy (we use the true value). 
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Figure 6. Top: Comparison between the observed CTI-free observations 
{{^fc}}i3.3 (black dots), and the modelled counts {{Aj;}}i3.3 (grey line) 
computed following the presented iterative procedure. The transits were 
generated from the typical reference image, at G = 13.3, with readout 
noise and background. 

Bottom: Residuals normalized by the noise: (A'^jfc — ^tk)/'^tk^ at the last 
stage of the image parameter estimation iterative procedure. 
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Figure 7. Normalized difference between the LSF model and the true under- 
lying flux distribution at different stages of the image parameter estimation 
iterative procedure (G = 13.3, readout noise and background, and typical 
reference image). The ideal case corresponds to a LSF model constructed 
from the observed data and the true image location. 



be equivalent to the quadratic sum of the photon noise and the read- 
out noise r. 

The agreement between observed and modelled counts 
(Fig. |6| does not significantly improve after 2 iterations, however 
the agreement between the LSF model and the reference image (see 
Fig.[7|, and the average image location bias as well as the location 
estimator standard errors (Sections |5.3| and \5.4\ can still improve 
after a certain number of iterations that essentially depends on the 
stellar magnitude. As a consequence, we stop the iterative proce- 
dure after a particular number of iterations that is determined for 
each magnitude beforehand. 



5 THEORETICAL AND ACTUAL LIMIT TO THE 
IMAGE LOCATION ACCURACY 

To be able to evaluate the accuracy of the image location esti- 
mation procedure, we first need to determine what the theoretical 
limit of any image location estimator is. This is done by comput- 
ing the Cramer-Rao bound (Section [5. 1^ , which shows that it de- 
pends uniquely on the image shape, flux, background and noise. 
Subsequently, and first in the absence of CTI, we verify that any po- 
tential bias of the Gaia image location estimator does not depend 
on the image location (Section \5.2\ . And then the estimator bias 
and standard errors are calculated as a function of G, image ref- 
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erence width, and for different operating conditions (Sections |5.3| 
and |5.4[ ). We compare the latter to the theoretical limit and eval- 
uate the efficiency of our estimation procedure in the absence of 
radiation damage. Then we estimate the irreversible loss of accu- 
racy intrinsic to radiation damage, by computing the Cramer-Rao 
bound for a 'damaged LSF' generated from the data set of damaged 
observations (Section [5. 5[ ). Ultimately, we apply the Gaia image 
location estimator to the damaged observations without any CTI 
mitigation. This allows us to characterize the radiation damage in- 
duced location bias (Section [5.6^ , and check the consistency of our 
CTI effects simulation by comparing our results to experimental 
test results (Section[5.7|(. 



^xlO"* - 



5.1 Definition of tiie astrometric Cramer-Rao bound 

For a dataset with a known underlying probability density func- 
tion the Cramer-Rao minimum variance bound theorem gives the 
minimum reachable variance of a free parameter using any estima- 
tion procedure. In the case of estimating the location k of a one- 
dimensional image containing Np detected photons, the Cr amer- 



Rao bound crjj can be expressed as follows jLindegren 



1978 



1 



\I'{x) 



-dx 



(12) 



I{x) + + r'')/Np 

with I{x) a normalized one-dimensional flux distribution of the 
image along x, (5 the background and r the CCD readout noise. 

5.2 Location independent error and standard deviation 

The iterative procedure on a transit (described in Section |4]( pro- 
vides us with an image location estimation for iteration n. One 
can compute the image location error 5^ by directly comparing «: 
to the true image location «:truc: 



5k = K — Ktr 



(13) 



Before averaging over all the image location estimates (or the cor- 
responding errors) of a particular magnitude (for a particular ref- 
erence image, window size, background level and readout noise 
value), one first needs to check that the error does not significantly 
fluctuate as a function of the relative location offset from the pixel 
grid. The latter is simply given by the collection of estimated im- 
age locations {k}g- Fig. [s] shows an example of this variation: 
each point corresponds to the average location error over 25 ad- 
jacent sub-pixel positions, the error bars represent the standard de- 
viation of the points with respect to this mean. At the last stage of 
the iterative procedure, the set of estimated locations {k*-"'} shows 
that there is virtually no significant systematic variation across the 
relative location offsets. A certain number of iterations (7 for the 
brightest and 2 for the faintest) is needed to remove the variation 
introduced during the procedure initialization by the Tukey's Bi- 
weight centroiding algorithm. 

Having established that there is no significant error as func- 
tion of the relative location offset from the pixel grid it is allowed 
to average over all the transits within a particular magnitude (for 
a particular reference image, window size, background level and 
readout noise value) to find the bias: 



1 



(14) 



with t going through all transits of the transit selection. We will 
indicate the average over all transits of a particular magnitude as 
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Figure 8. Image location error , as a function of the location offset 

from the pixel grid }q, for the brightest (top) and the faintest magni- 

tude (bottom) at the last stage of the iterative procedure (respectively 7 and 
2 iterations). 



{5k) a- For this subset of transits we can now also compute the 
corresponding standard deviation: 



T 



The statistical uncertainty of this standard deviation is: 



and the statistical uncertainty of the bias is: 

7f 



(15) 



(16) 



(17) 



Summarizing, we can for all transits of a particular magnitude 
quantify the location bias as: {Sk)g ± "(i^)' '^l^^ location stan- 
dard deviation, hereafter called location precision, as: a^.a ± ^cr^ ■ 
In the absence of any significant bias, the latter can be referred to 
as the location accuracy. 



5.3 CTI-free location bias results per magnitude 

Figure |9] shows the image location bias as a function of G for the 
three different reference image widths, a sky background level set 
to the average sky brightness and the Gaia CCD operating condi- 
tions regarding the readout noise value (4.35 e~) and the size of 
the telemetry windows in the along-scan direction (12 pixels for 
G <16 then 6 AL pixels G >16). This set of conditions, hereafter 
referred as to Gaia operating conditions, constitutes the most re- 
alistic case of our study and also the most unfavourable case for 
the image parameter estimation procedure. Yet, one can observe 
from Fig.|9]that the location bias, {Sk)g, for none of the magni- 
tudes exceeds the level of 5 milli-pixels (~ 0.3 mas). Moreover: 
within the uncertainty of our measurement, and for the three differ- 
ent image widths, (<5k)g does not significantly deviate from zero. 
Hence we can establish that the Gaia image location estimator is 
a bias-free estimator in the absence of radiation damage. Increas- 
ing the window size in the along-scan direction or reducing the 
readout noise has no significant effect on {5k)g- Only setting the 
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G 

Figure 9. Comparison of the location bias {5k) it t^(a^) for different ref- 
erence image widths as a function of G measured in the Gaia operating 
conditions (cf. Section [O) . The dashed line coiTesponds to a measure- 
ment realized in the same conditions for the typical image width but without 
background. The error bars correspond to the statistical uncertainty 

on the location bias (see Eq. |17| . Note that for readability a slight offset has 
been introduced on the G axis for the narrow and wide reference image re- 
sults. There are two ordinate axes, left is the location bias in units of pixels 
and the right in units of mas. The same holds for the following figures. 



background level to zero seems to slightly decrease the bias for the 
faintest magnitudes. This effect is illustrated by the dashed line in 
Fig.|9] 

5.4 CTI-free location accuracy per magnitude 

To evaluate the efficiency of our estimator, we compare the mea- 
sured standard errors, (Jk, to the astrometric Cramer-Rao bound, the 
theoretical limit to the image location accuracy of any bias-free es- 
timator (see Section [5T| . The comparison results are summarized 
in Table |3] for different values of G, image widths, window size, 
background levels, and values of CCD readout noise. In Fig. [TO] 
we compare the accuracy of the Gaia image location estimator in 
the Gaia operating conditions, the Cramer-Rao bound computed 
for the same reference image width and level of readout noise, and 
the requirements as presented in Table [T] As one can see, the Gaia 
image location estimator performs remarkably well. The estimator 
standard errors are always below the requirements and this for any 
reference image width. Also the standard errors, within the mea- 
surement statistical uncertainty vs^ , strictly follow the Cramer-Rao 
bounds at every signal level. Note how stringent the Gaia require- 
ments are: for the wide reference image, the actual and theoretical 
limits to the image location accuracy are very close to the required 
accuracy. 

As mentioned in Section [T| the targeted performance predic- 
tions of Gaia contain a margin of 20% to take into account unmod- 
elled on-ground calibration errors including for instance residual 
bias. In this context we consider an estimator efficient if its stan- 
dard errors are within 10% of the Cramer-Rao bound and thus not 
consuming more than half the margin. The Gaia estimator rigor- 
ously fulfills this criteria. This is illustrated in Fig.[TO](bottom) for 
the three reference image profiles: the ratio between the estimator 
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G 

Figure 10. Comparison between the actual and theoretical limit to the image 
location accuracy, as a function of G. Top: The continuous lines correspond 
to the Gaia image location estimator accuracy standard errors, '^o-^ , 

measured in the Gaia operating conditions for the three different reference 
images: naiTow (blue), typical (black) and wide (red). The associated er- 
ror bars corresponds to the statistical uncertainty on the location accuracy 
(see Eq. |16( . The dotted line represents the Cramer-Rao bounds computed 
for the same reference image widths (same colour coding), window size, 
background level, and readout noise value. Finally, the dashed and dotted 
(green) line shows the Gaia required image location accuracy. 
Bottom: The ratio between a 'he Cramer-Rao bounds and the asso- 
ciated error bars are depicted as a function G. For the three reference image 
widths, the relative deviation does not exceed 10%. The Gaia image loca- 
tion estimator can thus be considered efficient in the absence of radiation 
damage. 

standard errors and the Cramer-Rao bounds remain below 1.1 (i.e. 
10% relative deviation). As expected, in both the theoretical and 
actual cases, an increase in the image width is directly translated 
into a loss in location accuracy. This loss varies linearly with the 
image FWHM. Table[3]shows that increasing the readout noise, the 
background level, and/or decreasing the window size also increases 
the Cramer-Rao bound and the Gaia estimator standard errors. 

In the absence of radiation damage, we established that in re- 
alistic operating conditions and from bright to faint magnitudes, the 
Gaia image location estimator is bias-free, efficient, and performs 
within the requirements, with a high accuracy close to the theoreti- 
cal limit. It is now important to characterize in detail the impact of 
radiation damage on the image location uncertainty. 

5.5 Radiation damage intrinsic uncertainty increase 

Computing the Cramer-Rao limit (Eq. \12\ for a flux distribution 
including the CTI distortion and taking into account the charge 
loss allows to quantify this intrinsic uncertainty increase induced 
by the radiation damage. As one can observe from Fig. [2] the CTI 
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Figure 12. Image location bias, {&k)g,d, resulting from the radiation damage effects on tlie stellar image and the use of an image location estimator that 
does not take into account these effects. (5k )g d corresponds to the mean location en'or for all transits at a particular G. It is measured in the Gaia operating 
conditions (window size and readout noise value) for the three reference image widths and with no background (dashed line) and a background level set to 
the averaged sky brightness (continuous lines). The left figure shows the results obtained for a trap density of 4 traps pixel"-*^, representative of a CI delay of 
~27 s as demonstrated by the comparison to experimental results (crosses). The right figure shows the results obtained for a trap density of 1 trap pixel"^, 
representative of a CI delay of 1 s. The very small en'or bars correspond to the computed statistical uncertainty on the measured bias, jj^ (see Eq.jTT). The 
crosses represent the average relative location bias computed from experimental tests earned out during RC2, for a CI delay of ~27 s and 1 s and the lowest 
background level. The large experimental uncertainties are not shown, but as can be noticed the overall trend and amplitude of the measured bias is well 
reproduced by our model. 



induced distortion sharpens the image profiles and renders them 
more asymmetric but the charge loss significantly decreases the 
signal-to-noise ratio. The latter effect prevails and, at a given G, 
causes an increase in the image location uncertainty. To generate 
Ld, the damaged flux distribution, w e pro cee d in a similar fashion 



4.3 



and 



4.4 1. First we place 



to the construction of L (cf. Sections ■ 
each data point from the damaged observations at the right sub- 
pixel position to create an oversampled damaged profile. Then the 
over-sampled profile is fitted by the special quartic spline so that we 
can use an analytical representation. The resulting minimum vari- 
ances on the estimate of an image location acquired by a damaged 
CCD, and thus accounting for the CTI effects, are summarized in 
Table|4]for different image widths, background levels and levels of 
radiation damage. 

In the Gaia operating conditions, the relative intrinsic un- 
certainty increase (or accuracy loss) can be as large as 23% (see 
Fig.[TTJ for the highest trap density and 6% for the lowest. Here we 
recall that this drop in active trap density results from the use of a 
more frequent CI: from a CI period of 27 s to 1 s. In both cases, 
this increase is more pronounced for narrower stellar profiles and 
peaks at a signal level of G = 15.875. Then, due to the mitigat- 
ing effects of the SBC at lower signal levels, one clearly observes 
a flattening of the uncertainty increase. This illustrates the critical 
importance of the two hardware mitigation tools (see Section[3]4]l, 
which are the only mitigation countermeasures capable of reducing 
the CTI induced intrinsic loss of accuracy, by physically preventing 
the electron trapping and thus the image distortion and charge loss. 

The Cramer-Rao bound computed for the damaged flux dis- 
tribution now constitutes the maximum achievable accuracy by any 
unbiased image location estimator in the presence of radiation dam- 
age. Although the loss of accuracy can be quite large, Fig.|l 1 [shows 
that the Gaia requirements would still be fulfilled, if an image lo- 
cation estimator that is bias-free and efficient enough can be elab- 
orated (excluding the wide reference image and highest trap den- 



sity case). In the next section, in order to assess the efficiency of 
the Gaia image location estimator without CTI effects mitigation 
(Section|4]( in the presence of radiation damage, we directly apply 
it to the data set of damaged observations. 

5.6 Radiation induced image location bias 

In this section we are interested in exploring the consequences of 
not accounting for the CTI effects during the image location esti- 
mation. We thus apply the Gaia image location estimator as pre- 
sented in Section|4]to the data set of damaged transits. In this case 
the image distortion shall not be accounted for in the LSF model 
construction. This is achieved by using Lu, the LSF model gener- 
ated from the CTI-free transits. After applying the procedure, one 
eventually obtains an estimated location k o for each transit, which 
after subtraction of the true image location Ktruc, gives us the error 
Sk ■ After averaging for a particular magnitude and CCD operating 
conditions, we obtain the image location bias induced by the CTI 
effects as a function of signal level, (6,^)0, d- The location bias re- 
sults from the mismatch between the observed profile shape and the 
modelled LSF used to estimate the location. Only one iteration of 
the scheme from Fig. [5] is performed since the LSF model cannot 
be improved using the damaged counts without taking into account 
the CTI effects. 

The results for different image widths, window sizes, and 
background levels are summarized in Table [5] and depicted in 
Fig. 12 (left) for a trap density of 4 traps pixel" ^ and in Fig. 12 
(right) for 1 trap pixel" ^. The bias strongly varies as a function of 
G. In the Gaia operating conditions including background, the lo- 
cation bias reaches a maximum for C? = 15. The mitigation effects 
of the SBC is clearly noticed for 15 < G < 18 as the bias is either 
reduced or levels off. For G > 18, the background plays an impor- 
tant role in limiting the image distortion and reducing the bias as 
can be seen by comparing the dashed and solid lines. From these 
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Figure 11. Top: Comparison between the Cramer-Rao bounds computed 
from the original flux distribution (dotted lines) and the constructed dam- 
aged flux distribution for the two trap densities: 1 trap pixel^^ (dashed 
line), 4 traps pixel (continuous line). Two reference images are here con- 
sidered: the narrow (blue) and the wide (red). The window sizes, back- 
ground level and readout noise value correspond to the Gaia operating 
conditions. The computed required image location accuracy is also shown 
(dashed-dotted green line). 

Bottom: The relative intrinsic loss of accuracy induced by radiation dam- 
age as a function of G for the three reference images: narrow (blue), typical 
(black), wide (red), and the two trap densities: 1 trap pixel"^ (dashed line), 
4 traps pixel" ^ (continuous line). The relative loss of accuracy corresponds 
to the relative difference between the Cramer-Rao bound computed from 
the original flux distribution and the constructed damaged flux distribution. 
Note the important difference in loss amplitude, for the two different trap 
densities: a reduction in the active trap density (e.g., by the means of CI) is 
directly translated into a gain in location accuracy of a similar factor Sim- 
ilarly, the flattening of the intrinsic loss for G > 15 is due to the effect of 
the SBC. 



results we can conclude that in the presence of radiation damage, 
and without any attempt at any stage to correct or mitigate the CTI 
effects, the estimator is strongly biased. Indeed the image location 
can be shifted from a tenth of a pixel up to half a pixel for the fainter 
stars in the no-background case. When the estimator is applied to 
the damaged observations simulated with a background level set to 
the average sky brightness, the location bias is not as dramatic at 
low signal level. Nevertheless the image location bias for any im- 
age width and any signal level is constantly higher than ~0.1 pixels 
in the 4 traps pixel~^ case, which is not acceptable. Changing the 
telemetry window size has no significant effect on the location bias. 
As can be seen from Fig. [12] (right), for a shorter CI delay (or CI 



period), and thus less active traps (here 1 trap pixel" ^), the location 
bias is significantly lowered with a minimum level of ~0.02 pixels. 
It is interesting to note that decrease in bias is scaled by the same 
factor (~ 4) as the decrease in trap density. Regarding the required 
performance, for the faintest magnitude this level of bias might be 
acceptable in a limited amount of cases (e.g., the bluest stars). How- 
ever, in most cases, and especially for the bright stars this level of 
bias inevitably requires a software-based CTI mitigation scheme. 

5.7 Comparison with experimental data 

In order to check how representative the results obtained from syn- 
thetic data are in terms of the overall amplitude of the CTI effects 
and also fluctuation as a function of signal level. Fig. [T2] shows 
results obtained experimentally from RC2 ( |Georges|2 008 , Brown] 
]2009). In the experimental case the location bias does not corre- 
spond to an absolute image location bias since the true image loca- 
tion is by definition unknown. The presented bias is thus the relative 
location bias. It is computed by comparing the stellar transits over 
the irradiated part of the CCD and the same stellar transits over the 
non-irradiated part of the same CCD. Taking into account the dif- 
ferences between real and synthetic data, as well as experimental 
uncertainties, the overall agreement between the results obtained 
from the RC2 and our simulations is remarkable. The combined 
mitigating effects of the SBC and the background are also notice- 
able in the test data at low signal levels. Hence not only the am- 
plitude of the location bias for different CI delays (or densities of 
active traps) is reproduced by our model but also the overall bias 
evolution over a wide range of signal levels: 7 magnitudes. The 
simulations suggests that the illumination setup (and resulting PSF 
width) as well as slight differences in background light between 
experiments can have a significant impact on the measured CTI ef- 
fects. This may explain observed discrepancies between the results 
from different RCs and within a RC. 

5.8 Damaged location estimation standard errors 

Finally we show the resulting standard errors, (Jk^^ ^ function 
of G in Fig. [Ts] the standard errors are larger than the theoretical 
minimum variance, especially for intermediate magnitudes. For the 
most severe radiation level, the standard eiTors are larger than the 
requirements, and for the lowest radiation level the requirements 
are barely met; the mismatch between modelled and observed line 
spread function implies a broader spread in the locations estimated 
by the ML algorithm. This effect is less pronounced for the lowest 
level of radiation as the distortion, and thus the mismatch, is less 
important. The overall variance remains quite low as compared to 
the bias. Table l6] summarizes these results for the three different 
image widths. 

6 CTI EFFECTS MITIGATION 

Correcting for CTI is a complicated task and not only because the 
induced charge loss and distortion are considerable (Fig. ]2j. The 
trapping probabilities (e.g., [Prod'homme et ar]|201I | depend on 
the electron density and thus the CTI effects vary with G. This 
variation is not linear, in particular due to the presence of the SBC 
that mitigates the CTI effects only at low signal levels. This can be 
clearly observed from Fig. [12] from both simulations and experi- 
mental data. An important consequence is that the stellar core and 
wings (in the CCD serial direction) will not experience the same 
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Figure 13. Comparison between the Cramer-Rao bound in the presence of 
radiation damage (dotted hnes) and the standard en'ors (continuous Hnes) 
obtained by applying the Gaia image location estimation procedure to the 
set of damaged observations without any CTI mitigation (i.e. using the LSF 
model generated from the CTI-free observations). The Cramer-Rao bounds 
were computed and the standard eiTors measured for the typical reference 
image considering the Gaia operating conditions and two different levels of 
radiation damage: 1 trap pixel"^ (black) and 4 traps pixel"^ (red). In the 
presence of radiation damage, and without CTI effects mitigation, the Gaia 
image location estimator cannot be considered efficient anymore as its stan- 
dard errors deviate significantly from the Cramer-Rao bound (in addition to 
the estimator being biased c.f. Fig. |12^ . 



distortion. These different contributions to the global stellar im- 
age distortion will nevertheless be collapsed into a one-dimensional 
signal. In addition, the location bias and charge loss will not be re- 
peatable for a particular star or signal level as the CTI effects de- 
pend on the state (empty or filled) of the traps prior to the stellar 
transit. During the mission each star will transit on average ~72 
times over the focal plane of Gaia. For each of these stellar tran- 
sits the scanning direction of the satellite will differ, and thus also 
the CCD illumination history that determines the trap state. It is 
also likely that the trap density will have increased between two 
consecutive transits. The CI will play an important role here, not 
only by decreasing the active trap density but also by simplifying 
the illumination history by reseting it every 1 s. Finally, it is im- 
portant to note that, as already mentioned, Gaia's launch and first 
year of operations coincide with the predicted peak of the Sun's 
activity for the current solar cycle, and that none of the Gaia mea- 
surements will be free of radiation damage. This will strongly limit 
our knowledge of the exact instrument LSF/PSF in space. 



6.1 Potential alternative approaches 

Figure [T4| summarizes the Gaia data processing chain in three dif- 
ferent stages at which a different set of data is available: (i) the raw 
data, (ii) the intermediate data, (iii) the science data. Each set of 
data is further explained in the figure. Different ways of handling 
the CTI effects in this chain are possible, and the literature provides 
us with a handful of correction procedures for photometric, spec- 
troscopic, and (very rarely) astrometric measurements carried out 
in the optical or at X-Ray wavelengths. 
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Figure 14. From raw to science data: summary of the Gaia data processing 
chain in three stages (middle boxes). The top boxes contain the required 
models/assumptions to go from one stage to the other. The top arrows sym- 
bolize the feedback occurring at each stage which enables the iterative im- 
provements of both models and data. The bottom boxes give some details 
about the data delivered/used at each stage. In principle, each set of data 
can be empirically coiTected for CTI, however we choose not to perform 
any correction but to model the CTI distortion as part of the image parame- 
ter estimation procedure (see Section[6^and below). 



CTI correction at the level of the raw pixel data: one can cor- 
rect the raw pixel data to obtain artificial CTI-free data and perform 
the rest of the data processing using the corrected raw data. This 
constitutes one of the most common approaches, and has been suc- 
cessfully used to correct the CTI effects on HST data for instance. 
Its main advantage is that is minimizes the impact of CTI on the 
remaining data processing chain, the correction being performed 
very close to the source of the problem. Either the photo-electron 
count correction is directly performed by means of a parametric 
empirical or semi-empirical formula (e.g., [Goudfrooij & Kimble| 
|2002[ |Dolphin|2009| l that determines the CTI induced charge loss 
as a function of signal level, background, radiation dose, and source 
position on the CCD. Or it is performed by 'comparing' the dam- 
aged observation to a simulated observation, for which the dam- 
age is simulated by an empirical or physically-motivated analyti- 
cal forward model of the charge transfer and trapping (e.g. 
"tow 21)03 ; Mass ey et aL][20Tol [Anderson & Be din 2070| 



Bris- 



Bris- 



^tow ( 2003 1 provides a detailed comparison between direct empir- 
ical and model-based corrections: while the direct correction can 
only correct photometric and spectroscopic point source measure- 
ments, a model-based correction allows for astrometric correction 
of arbitrary complex sources (extended, binaries etc.). The latter 
is more complex, i.e. computationally intensive, but versatile and 
potentially more accurate. In principle, the model-based correction 
requires the generation of a synthetic undamaged observation to 
be subsequently distorted by the CTI model. However, as compre- 
hensively described in [Massey' s t al. (2010), and first proposed by 
[Bristow et al.[ ( |20()5] l, one can avoid this step and iteratively remove 
the CTI induced image distortion by subtracting actual and simu- 
lated observations, assuming, in a first step, that the actual damaged 
observation is the CTI-free input signal. This relies on the assump- 
tion that the CTI effects correspond to a slight perturbation around 
the true image. Although promising, the model-based correction 
of the raw data at the pixel level has only been tested against the 
empirical direct correction, and mostly for photometric and spec- 
troscopic data. Massey et aL|^2010| l go one step further and assess 
the astrometric correction induced shift as a function of signal level 
and distance from serial register. Although the correction performs 
as expected, the accuracy of such a method cannot be guaranteed 
yet due to the lack of reference or CTI-free data that prevents the 
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measurement of the method absolute bias and standard eiTors. On 
top of this uncertainty regarding the final accuracy of this method, 
two other considerations preclude the direct use of this approach 
in the Gaia data processing before more investigations. First, the 
noise properties of a corrected pixel value are no longer simple and 
may introduce hard-to-track effects in the image location estima- 
tion procedure, and subsequently in the astrometric global iterative 
solution (AGIS) that combines all observations to infer absolute 
astrometry for each observed object. In particular, the assumptions 
on which the maximum likelihood estimation of the image param- 
eters is based, namely that the individual samples are statistically 
independent and described by the Poissonian model (Section [4.2[ l, 
no longer hold for the corrected samples. Secondly, the lack of full 
frame data and the binning of most telemetry windows implies that 
we lack the information required to perform a full pixel-based cor- 
rection. 

CTI correction at the level of the intermediate data: at this 
level, the correction is performed thanks to a parametric ad-hoc 
model (e.g., [Rhodes et al.|2007[|Schrabback et al.|2010) . It offers 
the advantage of being simple and fast to apply, and once formu- 
lated the model should be relatively simple to calibrate. However, 
the elaboration of such model is not trivial. It first requires a care- 
ful study of the CTI effects on the parameters extracted from the 
raw measurements as a function of a finite number of pre-selected 
variables. Subsequent to this study, the dependency of the CTI in- 
duced bias on the pre-selected variables must be mathematically 
described for each estimated parameter of interest. It is not guar- 
anteed that such a mathematical formulation is possible and the re- 
sulting models have by definition no predictive power. In the case 
of Gaia, the CTI-induced image location bias and charge loss could 
be parametrized as function of the signal level, background, radi- 
ation dose (or observation time), source position on the CCD, and 
illumination history (or time since last CI). Fig. [T2] shows an ex- 
ample of the image location bias dependence on the signal level 
and background. Comparison between Fig. [T2] left and right, also 
provides additional information about the dependence on the time 
since last CI. Such an approach was studied by FADS Astrium, but 
does not constitute the current baseline approach of the Gaia Data 
Processing and Analysis Consortium (DPAC) as it cannot handle 
complex scenes but only single stars. 

CTI correction at the level of the science data: this last potential 
approach is the most impractical. It also requires a parametric ad- 
hoc model, most likely impossible to formulate as the CTI effects 
are too entangled at the level of the science data. Moreover, the cal- 
ibration of such approach would require the use of reference data, 
which in the case of Gaia will be mostly not available. 

6.2 A complete forward modelling approach 

Due to the complexity of the CTI effects and the extreme accu- 
racy required in the image location estimation, as well as for the 
reasons mentioned above, the DPAC adopted a forward modelling 
approach. Thus in contrast to the solution applied to HST data, no 
direct correction of the raw data shall be performed, essentially to 
preserve the simple noise properties and avoid arbitrary assump- 
tions. Instead, the true image parameters are estimated in an itera- 
tive scheme, in which each observation is ultimately compared to 
a modelled charge profile for which the distortion has been simu- 
lated through an analytical CTI model, a so-called charge distor- 
tion model (CDM). This approach is illustrated by the schematic 
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Figure 15. Forward modelling approach to CTI mitigation jLindegren] 
2008): a CTI-free sampled image is generated following the method ex- 
plained in Section |4] and is subsequently distorted by a fast analytical CTI 
effects model, so-called CDM (Charge Distortion Model). The distorted 
counts are then compared to the observed counts. In an iterative procedure, 
the scene, the LSF model, and the CDM parameters are successively im- 
proved. Note that if CDM takes as input a two-dimensional signal, the CDM 
output needs to be binned in the CCD serial direction before comparison to 
the observed counts. AGIS, the Gaia Astrometric Global Iterative Solution, 
uses the scene parameters to estimate the astrometric parameters. AGIS also 
provides an updated estimate of the scene parameters that is corrected for 
'nuisance' parameters such as the satellite attitude. 
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Figure 16. CDM parameter estimation procedure: a set of CDM input sig- 
nals is generated by sampling the current best LSF or PSF model using the 
best estimate of the scene parameters at a pailicular G. CDM simulates a 
set of damaged modelled observations that is subsequently compared to the 
set of observed counts corresponding to a particular G. A fitting algorithm 
provides us with a set of CDM parameters that optimizes the agreement 
between the CDM predictions and the observed counts. 



depicted in Fig. [15] where the modelled counts are now described 
as follows: 



Afe 



D [aL (A: - k) + /3 1 c, h] 



(18) 



with D the CTI distortion applied to the sampled image using 
CDM, c a set of CDM parameters (e.g., trap species characteristics, 
electron density distribution parameter), and h a set of parameters 
that describes the illumination history (the most obvious being the 
time since the last charge injection). 

As illustrated in the Fig. [15] the scene, the CDM, and the 
instrument (LSF/PSF) parameters are iteratively adjusted until 
the modelled counts D [{\k}] agree with the observation {A'^fc}. 
Fig. [T6] gives the details of the CDM parameter update. It is im- 
portant to note that the model LSF cannot be directly generated 
from the observations anymore as they are now affected by CTI. 
During the mission, the LSF model will thus be extracted from a 
LSF library composed partly by modelled LSFs and by a subset 
of observations: mostly the single bright stars that are the least af- 
fected by radiation damage (i.e. early mission data and/or observa- 
tions close to a charge injection). If the CDM and the instrument 
model are properly calibrated, the estimated scene parameters sub- 
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sequently used to determine the stellar parallaxes should be unbi- 
ased and free of CTI. Since no direct correction is performed the 
noise properties of the observation should remain dominated by the 
photon and readout noise and thus a complex contamination of the 
rest of the data processing chain and its products is avoided. This 
approach can handle arbitrarily complex scenes and offers the ad- 
vantage being in accordance with the general Gaia data processing 
principle of self-calibration. A similar approach was successfully 
used to handle CTI effects on photometric and spectroscopic X-ray 
measurements performed by Chandra jTownsley et al.|2000[|2002] 
[Grant et al.|2004p . 

In the following (Section [6.3[ l, we demonstrate the ability of 
the Gaia CTI mitigation approach to reach the best achievable im- 
age location estimation accuracy for damaged observations (Sec- 
tion |53|l in the case of an ideal CDM, and ideally calibrated LSF 
and CDM parameters. Then we assess the actual performance of 
this approach regarding the recovery the image location estimate 
bias (Section [6.6[ l and image flux estimate bias (Section [6.7[ l, using 
the current best CDM candidate jShort et al.|2010^ . 

6.3 Testing the forward modelling approach 

In a first step towards a more complete validation of our approach, 
we would like to ensure that this approach, if perfectly calibrated, 
enables an unbiased estimation of the image location with high 
enough precision. To do so we estimate the (unknown) scene pa- 
rameters for the set of damaged observations, in the case of an ideal 
CDM and ideally calibrated LSF and CDM parameters. This ideal 
case is simulated by using Ld, the damaged LSF (cf. Section 
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This is allowed because in this scheme the true LSF and CDM pa- 
rameters correspond to a model that is capable of fully explaining 
the image distortion and the charge loss in the damaged observa- 
tions. 

Figures [17] and [TS] show the location bias and the estimator 
standard errors obtained in these conditions, for the two different 
levels of radiation damage, and for the typical reference image. The 
results obtained for the two other reference images can be found in 
Tables[5]and[6] As one can see, in the case of ideal CTI mitigation, 
the location bias in the presence of radiation damage is now compa- 
rable to the one obtained for the CTI-free observations (see Fig.|9]l; 
the bias does not exceed 5 milli-pixels and does not significantly 
deviate from zero within the error bars (w(a^ ) , the statistical uncer- 
tainty), and this even for the most severe level of damage. Regard- 
ing the estimator standard errors, they comply with the Gaia re- 
quirements, even in the case of the most severe level of damage for 
most of the magnitudes. The bottom part of Fig. [18] shows that the 
location estimator including CTI mitigation performs efficiently for 
the lowest level of damage (i.e. less than 10% relative deviation). 
However it is interesting to note that even in this favourable case 
(ideally calibrated LSF model and CDM parameters), the relative 
deviation of the estimator precision from the best achievable one 
can reach 20% for the intermediate magnitudes and the strongest 
level of damage. 

From these results we can conclude that a forward modelling 
approach to CTI mitigation, as presented in the previous section, 
allows the recovery of the CTI-induced location bias and enables 
the bias-free estimation of the image location at the required pre- 
cision, close to the theoretical limit. This level of performance is 
achieved in the favorable conditions of a very good LSF model and 
the CDM parameter calibration, but for the strongest expected im- 
age distortion in the Gaia operating conditions, i.e. stars located the 
furthest away from the last CI and a density of traps equivalent to 
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Figure 17. Location bias {5k., d) i i'(i5^) after a CTI mitigation using the 
LSF model Lj^, representative of an ideally calibrated forward modelling 
approach. These results are obtained in the Gaia operating conditions, con- 
sidering the typical reference image only and for two different levels of 
radiation damage: 1 trap pixel" ^ (black) and 4 traps pixel" ^ (red). As one 
can see, in these conditions, the forward modelling approach allows for a 
full recovery of the CTI-induced location bias. Note that for readability a 
slight offset on the G axis as been introduced between the results for the 
two trap densities. 
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Figure 18. Comparison between the Cramer-Rao bound in the presence of 
radiation damage (dotted lines) and the measured standard errors (continu- 
ous lines) when applying the Gaia image location estimation procedure in- 
cluding CTI mitigation to the set of damaged observations. These results are 
obtained in the Gaia operating conditions, considering the typical reference 
image only and for two different levels of radiation damage: 1 trap pixel 
(black) and 4 traps pixel ~^ (red). The CTI mitigation is performed follow- 
ing the forward modelling approach, using Ljj (i.e. and ideally calibrated 
CDM and LSF model). In these conditions the estimator standard en'ors are 
below the requirements (green dashed and dotted line) even for the most se- 
vere damage (top). The relative deviation from the best achievable accuracy 
(bottom) is reasonable (below 1 0%) in the 1 trap pixel ~ ^ case but can reach 
20% in the 4 traps pixel^^ case for the intermediate magnitudes. 
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the predicted upper limit to the Gaia end-of-life accumulated radi- 
ation dose. 



6.4 Current best CDM candidate 

The elaboration and calibration of a CDM that allows to reach the 
level of performance presented in Section [63] is challenging. The 
presented mitigation scheme requires a CDM that must be both ac- 
curate and fast, as the iterative procedure is performed for each 
observation, and the CDM distortion applied at each iteration. The 
DPAC strategy regarding the elaboration of such CDM is detailed 
by |van Leeuwen & Lindegrenl^2007|> and|van Leeuwen|j2007^ and 
1 bylProd'homme ( 201 lb. In this study, to 



a short summary is given byl 
demonstrate the validity of our CTI mitigation approach including 
a CDM and thus assess its present actual performance, we use the 
current best CDM candidate (later referred as to CDM for simplic- 
ity) as it is described in [Short et al.| l (20T0} , and for which a first 
comparison of its outcomes to experimental test data is presented 
in |Prod'homme et al.|f2010) . 

CDM is based on the common Shockley Read Hall formalism 
( [Shockley & Read|1952"{|Hall|1952) and describes the capture and 
release processes in a statistical way. To cope with the computa- 
tional speed requirement, it suppresses the treatment of the numer- 
ous charge transfer steps required to transfer the signal from one 
CCD end to the other, but computes the signal transit in a single 
calculation making use of several assumptions jShort et al.|2010| l. 
CDM is able to simulate the CTI effects in TDI and imaging mode 
for any kind of signal (single, double stars, spectrum etc.). The 
CDM free parameters are: 7 which determines how the volume 
of the electron packet grows as electrons are added, /3 the back- 
ground light (respectively denoted /3 and Sdob in |Short et al.| ( |2010| , 
but changed herein for disambiguation), and three trap parameters 
per trap species, p, cr and r, respectively the trap density, the cap- 
ture cross-section, and the release time constant. It has to be noted 
that a more recent version of this model has been elaborated. This 
newer version incorporates a better handling of the charge injec- 
tion modelling and the possibility of simulating the serial CTI that 
occurs in the readout register. As charge injections are not explic- 
itly simulated in our synthetic data set and the serial CTI was not 
simulated, the hereafter demonstrated performances remain repre- 
sentative of the current performances of our mitigation scheme. 



6.5 The forward modelling approach initialization 

The iterative image parameter estimation procedure including CTI 
mitigation now involves three different sets of parameters to be suc- 
cessively improved: the scene, the LSF/PSF, and the CDM param- 
eters (see Fig.|15[l. Reaching a stable solution in these conditions is 
complex; each set of parameters needs to be initialized with values 
not too far off from the 'true' ones for the iterative procedure to 
converge. 

LSF model: as already mentioned, the LSF model cannot be gen- 
erated from the damaged observations, as they are not directly rep- 
resentative of the instrument anymore. During the mission the LSF 
model will partly be generated using the least damaged observa- 
tions of single bright stars. In the following we thus use Lu, the 
LSF model generated from the CTI-free observations at a particu- 
lar G. This constitutes a favourable yet realistic case. 



Scene parameters: the initial estimate for the image location k*-"' 
is determined using the Tukey's Biweight centroiding algorithm, 
the initial flux estimate a'"' corresponds to the sum of the observed 
counts after background subtraction (and as in the rest of the study 
the background /3 is considered to be known). Hence the initial lo- 
cation and flux estimates are biased by the CTI effects. However 
one should note that due to its construction the LSF model con- 
tains some information about the true location of the observations 
in its zero-point. This is still reasonable as we have so far ignored 
that during the mission the astrometric solution (AGIS) will pro- 
vide extra information about the true location of each observation 
through a feedback mechanism (see Fig.|15[l. 

CDM parameters: it is first important to realize that several fun- 
damental differences exist between CDM and the detailed Monte- 
Carlo CTI effects model that we used to simulate the damaged 
observations: the most important ones being related to the charge 
transfer simulation, the computation of the capture and release 
probabilities, and the modelling of the electron density distribu- 
tion. Hence, in this context, no 'true' CDM parameters exist but 
only CDM parameters that allow the reproduction of the simulated 
damaged observations. This actually constitutes a similar situation 
to the one that will be experienced during the operation of Gaia. In- 
deed, due to the simplifications intrinsic to the elaboration of a fast 
analytical model of a complex phenomenon, even with the right pa- 
rameterization, the agreement between the damaged observations 
and the CDM predictions will not be perfect. 

Furthermore: although the trapping occurs during the transfer 
of two-dimensional stellar images, only one-dimensional informa- 
tion is accessible from the binned observations. The CDM distor- 
tion can be applied to a one- or a two-dimensional CTI-free signal. 
In the latter case, one needs to reconstruct a PSF and the result- 
ing modelled counts must be binned prior to a comparison with 
the observed damaged counts. In our study we generally obtained 
a significantly better agreement between CDM predictions and the 
damaged observations by applying the CDM distortion to a one- 
dimensional signal. In the following we thus only present results 
obtained in this case. In reality this might be different, in partic- 
ular due to the serial CTI that was not taken into account here. 
Yet, if a comparable performance level can be achieved, the one- 
dimensional option would still be preferred during the mission for 
the ID binned data as it presents the advantage of saving a signif- 
icant number of computations. To obtain an initial set of CDM pa- 
rameters that describes reasonably well the damaged observations, 
we use Lu as input signal, and fit the CDM predictions to the dam- 
aged observations {{A^fe}}G for a particular G and set of operating 
conditions (i.e. windowing scheme, background and readout noise 

I between the 
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level). The fitting procedure minimizes the x (Eq. 
CDM predictions and the damaged observations. The fitted param- 
eters are 7, p, cr, and r (see Section [6!4) , /3 is fixed to the true value. 
At this stage, the fitting procedure is an evolutionary algorithnj^ 
that uses two mechanisms, mutation and cross-over. It is applied 
on an initial population of 100,000 parameter sets and evolves to- 
wards smaller generation after generation. After 10 generations, 
we select the set of parameters with the smallest . This set of pa- 
rameters can be further improved by using the downhill simplex 
minimization method l |Nelder & Mead|1965 ). Fig.[T9]gives an ex- 
ample of the obtained agreement between the CDM outcomes and 
the damaged 'observations' (generated with the Monte Carlo model 



http : / /watchmaker . uncommons . org/ 



Image location in the presence of radiation damage. 



17 




-6 -4 -2 2 4 6 



Along-scan coordinate [pixels] 

Figure 19. Top: Comparison between the CDM predictions after ini- 
tialization of the CDM parameters (grey) and the damaged observations 
{{^fe}}i3.3 (red dots) under the Gaia operating conditions and for a radi- 
ation level of 4 traps pixel" ^ . To enable the observation of the CTI induced 
distortion the input signal Lu is also shown (black line). 
Bottom: Residuals normalized by the photon noise, the reduced is ~3.0. 

described in Section[3| at a particular value of G. This example is 
representative of the best level of agreement achieved after apply- 
ing the described initialization procedure. The illumination history 
parameters, h, will be fixed to the reconstructed illumination his- 
tory. Here h is only the time since last CI that is set to infinity as 
no CI has been explicitly simulated. The effect of not calibrating 
for disturbing stars, i.e. stars located between the last CI and the 
star of interest, will be studied in the second part of this study. It 
can however already be mentioned that stars located between a CI 
and the star of interest are only disturbing if they are located in the 
same pixel column (or an adjacent one) and that, for a CI period of 
1 s, the number of disturbing stars is expected to be very low even 
for the densest parts of the sky ( |Holl et al.|20Ilb^ . 

6.6 Image location bias and accuracy recovery 

As mentioned in the previous section, after initialization, the CDM 
parameters can still be further improved by the use of the downhill 
simplex method. In the following we thus distinguish between two 
different cases: {a) the CDM parameters have been fully optimized 
and the scene parameters are then estimated, no more iterations are 
performed, {b) the CDM parameters have not yet been fully opti- 
mized and are refined as part of the image location iterative pro- 
cedure (Fig. |I5[ ). In the latter case, once each set of parameters is 
initialized, the iterative procedure is performed as follow: (i) the 
scene parameters are first estimated using Lu and the CDM pa- 
rameters 7'°', p'°\ o-'°\ r*°^ and P, then (ii) the CDM parameters 
are updated as presented in Fig.[T6]using the newly estimated scene 
parameters, and (iii) a new scene parameter update is performed. In 
our study, the CDM parameter update (Fig.|I6[> is performed using 
the downhill simplex minimization method (N elder & Mead|I965^ 
only; as we shall see it proves to be quite inefficient at this stage of 
the procedure. A maximum likelihood based procedure would be 
better suited and is currently being developed to perform this task 
in the Gaia data processing. 
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Figure 20. Image location bias recovery using the CTI mitigation forward 
modelling approach including the current best CDM candidate. These re- 
sults are obtained in the Gaia operating conditions, considering the typ- 
ical reference image only, for two different levels of radiation damage 
(1 trap pixel"^ in black, and 4 traps pixel"^ in red), and for two different 
initial optimizations of the CDM parameters: fully optimized (continuous 
line) and iteratively improved (dashed line). For comparison the average 
location bias measured for the CTI-free observations is also shown (blue 
dotted line). The CDM parameters are calibrated per magnitude and as a 
result there is a different set of CDM parameters for each magnitude and 
the resulting final agreement between modelled and observed counts varies 
from one signal level to the other This explains the bias oscillations (in 
particular for the highest trap density). 




13 14 15 16 17 18 19 20 



Figure 21. Image location standard errors resulting from the use of the CTI 
mitigation forward modelling approach including the current best CDM 
candidate for two different levels of radiation damage (1 trap pixel ~^ in 
black, and 4 traps pixel^^ in red), and for two different initial optimization 
of the CDM parameters (fully optimized as a continuous line, and iteratively 
improved as a dashed line). The lower panel shows the relative deviation 
from the best achievable accuracy. The Cramer-Rao bound in the presence 
of radiation damage and the Gaia requirements are also shown (respectively 
the dotted lines and the preen d^i.^hed and Hotted line^. 
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Figure |20] shows the remaining image location bias after us- 
ing the CTI mitigation forward modelling approach including the 
current best CDM candidate for the two different initial optimiza- 
tions of the CDM parameters. These results should be compared 
to Fig. [12] that shows the location bias when no CTI mitigation 
is applied, and to Fig. [TT] that shows the ideal performance of the 
presented mitigation scheme. All these results are summarized in 
Table |5] for the three different reference images. The current best 
CDM candidate does not allow for a total recovery of the location 
bias at each magnitude, however this bias is considerably reduced 
for both levels of radiation damage. For instance, in the case of 
the lowest trap density and an optimal CDM optimization, the bias 
does not exceed the level of 0.005 pixels, while without any miti- 
gation and in the same conditions the bias reaches 0.05 pixels (see 
Fig.[T2j. Fig.|20|shows that at the faint end the bias is significantly 
lowered by a better optimization of the CDM parameters, and that 
the CDM parameter update performed during the iterative proce- 
dure in these conditions is too limited to recover the full potential 
of CDM. This means that the CDM initialization is a crucial step 
of the CTI mitigation scheme, especially at faint magnitudes, and 
if the CDM parameters are iteratively refined, the simplex method 
seems not to be efficient enough. The bias oscillations as a func- 
tion of G are due to the fact that the CDM parameters have been 
calibrated per magnitude. As a result a different level of agreement 
between observed and modelled counts is achieved for each signal 
level. 

Figure [27]shows the measured location standard errors for the 
typical reference image (see Table |6] for the other reference im- 
ages). As expected (see Fig. [T8]( for the most severe level of ra- 
diation damage, the Gaia requirements are not met for the bright 
magnitudes, and the relative deviation from the best achievable ac- 
curacy is large: it almost reaches 100% at G = 14.15. However 
when considering the more realistic trap density of 1 trap pixel" ^, 
the standard errors are safely below the requirements, and the rel- 
ative deviation from the Cramer-Rao bound remains below 10%. 
In these conditions our estimator thus remains efficient, even if, as 
already explained, biased. Finally, it has to be noted that the overall 
precision of the location estimation seems to be quite insensitive to 
the fine tuning of the CDM initial parameters. 



6.7 Image flux bias recovery 

This study focused on the estimation of the image location param- 
eter as it is the most critical image parameter to be determined 
for astrometry. However the accurate estimation of the integrated 
image counts is also important as this forms the basis for the G- 
band photometry. The photometry, when combined with the par- 
allax measurements, will provide the absolute luminosities of the 
stars observed by Gaia. In addition to this fundamental parame- 
ter the multiple observations of each source constitute an all sky 
variability survey, providing another treasure trove of astrophysi- 
cal information. For both applications high photometric accuracy, 
complementary to the astrometric accuracy is required. We refer 
to |Jordi etal.1 (201Ol for more details on Gaia's photometric capa- 
bilities. As can be observed from Fig. |2] CTI induces not only an 
image distortion but also an important charge loss that, if not prop- 
erly taken into account, biases the image flux estimation. In order 
to judge the capability of our approach to CTI mitigation to achieve 
high photometric accuracy we show in Fig.|22]the image flux esti- 
mation bias (in units of magnitude) achieved with and without CTI 
mitigation using CDM. These biases are compared to the photo- 
metric performance predictions in |Jordi et al.| ( |2010l l including the 
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Figure 22. Comparison between the image flux estimation bias (in units 
of magnitude) before (top) and after CTI mitigation (bottom) in the Gaia 
operating conditions considering the typical reference image for two differ- 
ent levels of trap density: 1 trap pixel" ^ (black) and 4 traps pixel" ^ (red). 
These results are compared to the Gaia requirements fJ ordi et al.|201 0't for a 
transit across a single CCD (green dashed and dotted line). In the case of the 
CTI mitigation (bottom) we considered two different initial optimizations 
of the CDM parameters: fully optimized (continuous lines) and iteratively 
improved (dashed lines). 



intrinsic loss in photometric precision induced by CTI ( |Jordi et aT] 
|20I0[ l, where the numbers in their figure 19 have been translated to 
the photometric errors expected for a transit across a single CCD. 
The 'safety margin' in equation (6) in |Jordi et al.| ( |20I0) was omit- 
ted. 

Despite the rather strong bias induced by CTI in the image flux 
estimation (see Fig.|22]top), the presented CTI mitigation scheme 
allows to eliminate most of it (see Fig.|22]bottom). It thus allows an 
unbiased estimation of the image flux within the requirements if the 
estimation procedure is precise enough. In this context the current 
CDM performances are remarkable: for the highest trap density and 
when no mitigation is applied, the flux bias can reach ~0.25 magni- 
tudes at G = 15.875, after mitigation we measure a flux estimation 
bias of 0.0029 magnitudes for the same magnitude, well below the 
requirement of 0.0067 magnitudes. It is also interesting to note that 
the image flux estimation is much less sensitive to the calibration of 
CDM than the image location estimation. Indeed Fig. |22] (bottom) 
shows that a similar level of performance is obtained for a fully 
optimized or an iteratively improved CDM calibration. 



7 DISCUSSION 

Throughout the paper we assumed a fully functional SBC. Never- 
theless the manufacturing of a SBC is a complex process and the 
CCDs of the Wield Field Camera channel of the Advanced Cam- 
era for Surveys on board HST have been reported not to contain 
the SBC present in their design ( [Anderson & Bedin|20"T0} . Fig. [23] 
shows the importance of the SBC CTI mitigation for achieving the 
Gaia requirements at low signal levels. If the SBC was not present 
(dashed line), the intrinsic loss of accuracy at low signal levels 
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Figure 23. Top: The relative intrinsic loss of accuracy induced by radiation 
damage as a function of G and as computed in Section [53] The continuous 
line is the same as the black dashed line in Fig. |ll| and corresponds to the 
case where the SBC is functional. The dashed line is obtained using simu- 
lated transits for a CCD containing no SBC. Note the difference in ordinate 
scale range with Fig. | 11 [ (bottom). 

Bottom: Absolute image location bias as a function of G and as computed 
in Section l5!6l The continuous line is the same as the black continuous line 
in Fig. ^[(typical image width, Gaia operating conditions, 1 trap pixel^^) 
and corresponds to a functional SBC case. The dashed line was obtained for 
simulating star transits in the same conditions but with a CCD containing 
no SBC. 



would reach 200% instead of only 5% in the presence of a SBC 
(continuous line and Fig.[TT] bottom). The extra location bias in- 
duced by a missing SBC is also significant although there is no 
particular reason for which the presented CTI mitigation scheme 
would not be able to recover it. In the case of Gaia the SBC has 
been demonstrated to be functional using experimental tests. For 
instance, Fig. 5 in |Prod'homme et ar] l |201iy shows the effect of the 
SBC on the CTI-induced fractional charge loss as a function of sig- 
nal level. However a recent study by ^Kohley et al.|p009^ identified 
a non-functional SBC in the upper half of a Gaia CCD. Based on 
more tested devices, Seabroke et al. (in prep.) will show that a sig- 
nificant number of the Gaia CCDs could be affected by this issue. 
Using the same methods as in this paper, they will evaluate the ex- 
tra loss of accuracy and location bias induced by a non-functional 
SBC in the CCD upper half only; in contrast to the no SBC case in- 
vestigated here, Seabroke et al. will show that there is only at most 
10% extra loss of accuracy. 



8 CONCLUSIONS 



curacy. The underlying principle of this study consisted of a sys- 
tematic comparison between the computed theoretical limit to the 
image location accuracy, the Cramer-Rao bound, and the actual per- 
formance of the Gaia image parameter estimation procedure under 
realistic Gaia operating conditions. The image location estimation 
bias and the associated standard errors were measured by applying 
the Gaia image parameter estimation procedure to a large set of 
synthetic data accounting for different stellar image widths, mag- 
nitudes, and background levels. We considered two different ac- 
tive trap densities; they are representative of two different levels 
of hardware mitigation in the presence of a radiation dose equiv- 
alent to the upper limit to the expected end-of-mission accumu- 
lated dose. The lowest active trap density constitutes the most re- 
alistic case because of a lower expected radiation dose and more 
importantly because of the planned injection of artificial charges in 
the CCD every second. In this context, a total of 41,472 synthetic 
two-dimensional Gaia-like observations have been generated using 
a detailed Monte-Carlo model of the CCD charge collection and 
transfer, and the radiation induced trapping. The dataset is readily 
available for the Gaia scientists to continue to test and further im- 
prove these critical steps in the Gaia data processing that are the 
image parameter estimation and the mitigation of the CTI effects. 
The main conclusions we can draw from this study are: 

The Gaia image location estimation procedure is bias-free and 
efficient in the absence of radiation damage. We showed that 
under realistic operating conditions, and from bright to faint mag- 
nitudes, the Gaia location estimator performs within the require- 
ments at an accuracy close to the theoretical limit for CTI-free ob- 
servations. 

The radiation damage effects induce an irreversible loss of ac- 
curacy that is independent of any image location estimator. It 

can only be avoided by the use of hardware CTI countermeasures 
that physically prevent the trapping. In the theoretical limit (i.e. 
perfect CTI calibration at the image processing level), the loca- 
tion accuracy loss is still acceptable when compared to the Gaia 
requirements: it can reach 6% in the lowest trap density case and 
24% in the highest case. Due to the presence of a supplementary 
buried channel in each of the Gaia CCD pixels, the accuracy loss 
stops increasing for stars fainter than G ^ 16. 

A CTI mitigation procedure is critical to achieve the Gaia re- 
quirements. We showed that if CTI is not taken into account in 
the image parameter estimation procedure, the resulting image lo- 
cation estimations are significantly biased. In the Gaia operating 
conditions, the most important bias is obtained for the widest type 
of stellar images at G = 15: 0.05 pixels in the lowest trap density 
case and 0.2 pixels in the highest case. For comparison, at this mag- 
nitude the requirement on the image location accuracy for a G2V 
type star is 0.0045 pixels, at least an order of magnitude smaller 
than the measured bias. 

The CTI-induced image location bias varies significantly with 
the stellar image width and the background level. This is par- 
ticularly relevant for experimental studies in which the image flux 
distribution and the background level cannot be absolutely known. 
At faint magnitudes small differences in the experimental setup can 
lead to significant differences in the measured CTI effects. 



We have presented a detailed characterization and evaluation of the 
impact of CCD radiation damage on the Gaia image location ac- 



In principle, a complete forward modelling approach to CTI 
mitigation allows for an accurate and bias-free estimation of the 
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true image location from a damaged observation. We demon- 
strated that the forward modeUing of a damaged observation using 
ideally calibrated models for both the PSF/LSF and the CTI ef- 
fects provides a location estimate that on average never exceeds 
0.003 pixels and does not deviate from zero within the error bars 
from bright to faint stars and for the two considered levels of ra- 
diation damage. The accuracy reached using this CTI mitigation 
scheme complies with the Gaia requirements. In the case of the 
lowest active trap density, this method even allows for the recov- 
ery of the theoretical limit to the image location accuracy in the 
presence of radiation damage. 

H calibrated well enough, the current best candidate for the 
charge distortion model (CDM) associated with the forward 
modelling approach allows significant image location and flux 
estimate bias recovery. In these favourable conditions (simple 
illumination history, 1 trap species, no serial CTI, well calibrated 
LSF model, and close to optimal CDM parameters) yet for a trap 
density level representative of the end-of-mission accumulated 
radiation dose, the Gaia image location accuracy is preserved. 
In the Gaia operating conditions and after CTI mitigation using 
the current best CDM at our disposal, the maximum measured 
location bias is ~0.005 pixels for the lowest radiation level and 
~0.017 pixels for the highest. 



as function of magnitude, shape of the stellar profile, illumination 
history (time since last charge injection), and mission time (or trap 
density). These models will be used to disturb the image locations 
(observation times) processed by AGIS. The study of the resulting 
astrometric parameters will then allow us to characterize and eval- 
uate the impact of CCD radiation damage on Gaw's astrometry. 
Also the effect of disturbing stars between a charge injection and a 
target star will be assesed. 

The future steps outlined are crucial ingredients in the suc- 
cessful radiation damage mitigation strategy for Gaia, enabling the 
extraction of the best scientific performance from this exciting and 
much anticipated mission. 
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9 FUTURE WORK 

Estimating the location of an image to milli-pixel accuracy is an ex- 
tremely challenging exercise, in which no detail must be neglected. 
This is especially true in the presence of radiation damage as shown 
in this paper. The work presented here is not the final word on the 
Gaia image parameter estimation procedure. Indeed the estimation 
procedure must be tested and improved further using synthetic and 
experimental data. In addition the elaboration and the calibration 
of the charge distortion model is a key element in the success of 
the presented CTI mitigation scheme. In this study, we have estab- 
lished the level of agreement with the damaged observation that any 
CDM must achieve to recover a bias-free image location estima- 
tion. We intend to test if the current best CDM candidate is capable 
of reaching such agreement with experimental data and improve it 
if necessary. Regarding the calibration of the radiation damage pa- 
rameters, the periodic charge injections will enable us to monitor 
and characterize the radiation damage during the mission. In addi- 
tion the charge injections will act to reset the illumination history. 
We intend to study what are the parameters that one can infer from 
the study of the CTI effects on the charge injection signal, and how 
one can use these parameters to initialize and calibrate the charge 
distortion model. 

This study constitutes the first step in evaluating the impact 
of the CCD radiation damage on the final astrometric accuracy of 
Gaia. It is indeed not yet clear in detail how a biased and less pre- 
cise estimation of the image location, as induced by CTI, prop- 
agates into the astrometric parameters derived by the Gaia as- 
trometric global iterative solution (AGIS). In the follow-up paper 
( |HoII et ar]|2011a] l we investigate this particular question by us- 
ing a small-scale version of AGIS, AGISLab ( |HolI et al.|[2509l l, 
that will allow us to perform a careful error propagation analysis 
for different cases (no CTI mitigation and optimal mitigation ver- 
sus CTI-free). The study presented here will be used to construct 
for each case a model that provides location bias and uncertainty 
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2. List of acronyms used in this paper. 



Acronym 


Definition 


AC 




AF 






rtftLlUlllclllC VJHJUcU lUCIallVC ijUiUllUIl 


AL 


ALong scQn 


CCD 


X^llCligC ^^-ULIUICU. L'CVlV^C 


CDM 


1^ n nTfTf* Til cfjTrfi jtti ixAr^ri^*! 


CFMnA 

W JJ/IVIVJ /V 


PTT FffppU MnHpk frtr fJAia 


CI 


r^harop Tniprtinn 


CTI 


Charge Transfer Inefficiency 


ELSA 


European Leadership in Space Astrometry 


FWHM 


Full Width Half Maximum 


HST 


Hubble Space Telescope 


L2 


Lagrangian Point 2 


LSF 


Line Spread Function 


mas 


milli-arcsecond 


/las 


micro-arcsecond 


ML 


Maximum-Likelihood 


PSF 


Point Spread Function 


RC 


Radiation Campaign 


SBC 


Supplementary Buried Channel 


TDI 


Time-Delayed Integration 
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Reference image type 


typical 






typical 




typical 




typical 


Background level (e^pixel^-^ s~^) 














0.44698 




0.44698 


Readout noise {e~ ) 









4.35 




4.35 




4.35 


Window size (pixels) 


40 






40 




40 




18 


Trap density (traps pixel" ^) 






















Magnitude 


Cramer-Rao bound 
















(G-band) 


(10"'' pixel) 


(10 pixel) 














13.3 


1.67 


1.668 ± 0.075 


1.67 


1.672 ± 0.075 


1.68 


1.655 ± 0.074 


1.67 


1.662 ±0.074 


14.15 


2.48 


2.49 ± 0.1 1 


2.48 


2.50 ± 0.1 1 


2.48 


2.46 ±0.11 


2.48 


2.45 ±0.1 1 


15.0 


3.66 


3.69 ±0.17 


3.67 


3.71 ± 0.17 


3.67 


3.77 ± 0.17 


3.67 


3.75 ±0.17 


15.875 


5.48 


5.95 ± 0.27 


5.50 


6.00 ± 0.27 


5.50 


5.45 ± 0.24 


5.51 


5.63 ±0.25 


16.75 


8.29 


7.65 ± 0.34 


8.33 


7.65 ± 0.34 


8.34 


8.46 ± 0.38 


8.30 


8.34 ±0.37 


17.625 


12.40 


12.65 ± 0.57 


12.54 


12.73 ± 0.57 


12.57 


12.63 ± 0.56 


12.61 


12.70 ±0.57 


18.5 


18.55 


17.88 ± 0.80 


18.99 


18.03 ± 0.81 


19.09 


18.84 ± 0.84 


19.37 


19.34 ±0.86 


19.25 


26.29 


25.87 ± 1.16 


27.49 


26.56 ±1.19 


27.62 


27.83 ± 1.24 


28.53 


28.86 ±1.29 


20.0 


37.06 


37.12 ± 1.66 


40.26 


39.99 ± 1.79 


40.99 


38.08 ± 1.70 


43.40 


39.08 ±1.75 


Reference image type 


typical 






narrow 




typical 




wide 


Background level (e^pixel^-^ s^^) 









0.44698 




0.44698 




0.44698 


Readout noise (e^ ) 









4.35 




4.35 




4.35 


Window size (pixels) 


telemetry* 




telemetry 




telemetry 




telemetry 


Trap density (traps pixel" ^) 






















Magnitude 


Cramer-Rao bound 
















(G-band) 


(10^3 pixel) 


(10^3 pixel) 














13.3 


1.67 


1.680 ±0.075 


1.50 


1.416 ±0.063 


1.67 


1.658 ± 0.074 


1.95 


1.920 ±0.086 


14.15 


2.48 


2.54 ±0.11 


2.23 


2.25 ±0.10 


2.48 


2.46 ±0.11 


2.88 


2.83 ±0.13 


15.0 


3.66 


3.69 ±0.17 


3.30 


3.52 ±0.16 


3.67 


3.74 ±0.17 


4.27 


4.52 ± 0.20 


15.875 


5.48 


5.98 ± 0.27 


4.95 


5.47 ± 0.24 


5.51 


5.59 ± 0.25 


6.42 


5.69 ± 0.25 


16.75 


8.29 


7.69 ± 0.34 


7.50 


7.55 ± 0.34 


8.38 


8.45 ± 0.38 


9.85 


9.80 ± 0.44 


17.625 


12.40 


12.66 ± 0.57 


11.39 


11.13 ±0.50 


12.70 


12.94 ± 0.58 


14.94 


15.55 ± 0.70 


18.5 


18.55 


17.99 ±0.80 


17.44 


16.88 ±0.76 


19.49 


19.39 ±0.87 


22.99 


25.09 ± 1.12 


19.25 


26.29 


25.87 ± 1.16 


25.48 


26.33 ± 1.18 


28.65 


29.53 ± 1.32 


33.73 


31.20 ± 1.40 


20.0 


37.06 


37.12 ± 1.66 


38.43 


39.73 ± 1.78 


43.57 


41.66 ± 1.86 


51.11 


51.94 ±2.32 



Table 3. Comparison between the theoretical and the actual limit to the image location accuracy in the absence of radiation damage as a function of G, sky 
background, readout noise, image width, and size of the telemetry windows in the along-scan direction. The theoretical limit corresponds to the Cramer-Rao 
bound, and the actual to the Gaia image location estimator standard errors with v^^ the statistical uncertainty (Eg. |16| . While the image width has a 
significant impact on those limits (e.g., for a 20% increase in FWHM, one can note a ~25% decrease in accuracy for the faintest magnitude), the window size, 
readout noise, and background level only slightly affects the image location accuracy. 

* 'telemetry size' refers to the size of the windows as they will be transmitted to the ground segment during the operational phase of Gaia ^Paulet|2009^ : 
12 pixels in the along-scan direction for G < 16, and 6 pixels for G > 16. 



Reference image type 




typical 






naiTow 






typical 






wide 




Background level (e^pixel^^ s 











0.44698 






0.44698 






0.44698 




Readout noise (e") 




4.35 






4.35 






4.35 






4.35 




Window size (pixels) 




telemetry 






telemetry 






telemetry 






telemetry 




Trap density (traps pixel" ^) 





1 


4 





1 


4 





1 


4 





1 


4 


Magnitude (G-band) 










Cramer-Rao bound (10 


pixel) 










13.3 


1.67 


1.69 


1.72 


1.51 


1.52 


1.55 


1.68 


1.69 


1.73 


1.95 


1.96 


2.01 


14.15 


2.48 


2.53 


2.69 


2.23 


2.27 


2.40 


2.48 


2.53 


2.70 


2.88 


2.94 


3.15 


15.0 


3.67 


3.81 


4.28 


3.30 


3.43 


3.87 


3.67 


3.82 


4.31 


4.27 


4.45 


5.02 


15.875 


5.50 


5.79 


6.61 


4.95 


5.25 


6.10 


5.51 


5.82 


6.67 


6.42 


6.73 


7.64 


16.75 


8.33 


8.72 


9.97 


7.53 


7.98 


9.23 


8.41 


8.85 


10.06 


9.97 


10.43 


11.68 


17.625 


12.54 


13.08 


14.89 


11.44 


12.01 


13.69 


12.74 


13.35 


15.01 


15.11 


15.70 


17.46 


18.5 


18.99 


19.71 


22.16 


17.50 


18.22 


20.51 


19.54 


20.28 


22.59 


23.24 


23.98 


26.47 


19.25 


27.49 


28.54 


33.12 


25.54 


26.78 


30.51 


28.72 


29.93 


33.84 


34.01 


35.50 


39.57 


20.0 


40.26 


43.27 


53.28 


38.51 


40.27 


45.96 


43.64 


45.43 


51.19 


51.43 


53.54 


60.32 



Table 4. Comparison between the Cramer-Rao bounds computed for different image widths, background levels, and radiation damage levels. This comparison 
allows to characterize the intrinsic loss of precision and ultimately accuracy induced by radiation damage. This loss is relatively more important for the 
narrowest image, and increase with the trap density as expected. 
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Reference image type narrow typical wide 



Bacltground level (e^pixel^^ s^^) 


0.44698 


0.44698 


0.44698 


Readout noise (e^) 


4.35 


4.35 


4.35 


Window size (pixels) 


telemetry 


telemetry 


telemetry 


Trap density (traps pixel^^) 


1 


4 


1 


4 


1 


4 


Mitigation 


none 


none 


none 


Magnitude (G-band) 




Location bias {5k.)g 


±v^s^-^ (10-3 


pixel) 




13.3 


24.85 ±0.11 


93.48 ± 0.27 


29.06 ±0.11 


110.84 ±0.22 


35.64 ±0.13 


137.18 ±0.19 


14.15 


35.79 ±0.14 


137.10 ±0.22 


40.93 ±0.17 


156.83 ± 0.29 


49.68 ± 0.20 


188.40 ±0.28 


15.0 


43.54 ± 0.24 


166.65 ± 0.66 


48.53 ± 0.23 


180.64 ±0.47 


56.48 ± 0.27 


207.34 ± 0.54 


15.875 


44.06 ± 0.44 


161.29± 1.11 


46.81 ± 0.43 


170.11 ±0.82 


51.19 ±0.48 


187.06 ± 0.74 


16.75 


37.70 ± 0.56 


139.54 ± 1.12 


39.75 ± 0.60 


143.66 ± 0.93 


43.03 ± 0.64 


157.81 ± 0.97 


17.625 


31.26 ±0.81 


120.41 ± 1.18 


33.02 ± 0.87 


126.87 ± 1.04 


33.66 ± 0.99 


139.18 ± 1.20 


18.5 


26.09 ± 1.07 


108.10 ± 1.29 


26.11 ± 1.28 


114.60 ± 1.48 


33.79 ± 1.58 


130.81 ± 1.50 


19.25 


28.67 ± 1.64 


105.25 ± 1.92 


31.02 ± 1.81 


109.66 ± 2.03 


25.92 ± 2.14 


119.33 ±2.20 


20.0 


24.92 ±2.51 


84.61 ± 2.73 


19.97 ± 2.86 


87.34 ±3.15 


24.16 ± 3.57 


95.08 ± 3.71 


Mitigation 


ideal 


ideal 


ideal 


Magnitude (G-band) 




Location bias (Sk)g 


±V^s^-^ (10-3 


pixel) 




13.3 


0.01 ±0.10 


0.02 ±0.12 


-0.05 ±0.11 


0.18 ±0.12 


0.09 ±0.13 


-0.19 ±0.13 


14.15 


-0.16 ±0.14 


0.05 ±0.16 


0.15 ±0.16 


0.11 ±0.20 


-0.17 ±0.20 


-0.15 ± 0.20 


15.0 


0.27 ± 0.22 


0.55 ± 0.36 


-0.32 ± 0.25 


0.29 ± 0.33 


-0.25 ± 0.28 


0.11 ±0.32 


15.875 


-0.12 ± 0.36 


-0.30 ± 0.52 


0.12 ±0.38 


1.03 ±0.47 


-0.33 ± 0.45 


0.44 ± 0.47 


16.75 


0.00 ± 0.50 


0.51 ±0.63 


-0.14 ±0.58 


0.40 ± 0.62 


-0.07 ± 0.65 


0.81 ± 0.70 


17.625 


-1.63 ±0.77 


0.25 ± 0.87 


0.22 ± 0.85 


0.00 ± 0.90 


-1.82 ±0.97 


1.04 ± 1.10 


18.5 


-0.54 ± 1.07 


-1.94 ± 1.24 


0.90 ± 1.26 


-1.52 ± 1.42 


2.04 ± 1.58 


-0.38 ± 1.51 


19.25 


-0.28 ± 1.66 


-1.98 ± 1.91 


-2.09 ± 1.86 


-0.67 ± 1.95 


1.30 ±2.15 


-0.75 ± 2.18 


20.0 


4.84 ± 2.44 


0.30 ±2.91 


-3.08 ± 2.96 


-2.06 ±3.18 


-2.63 ± 3.53 


6.64 ± 3.58 


Mitigation 


CDM 


CDM 


CDM 


Magnitude (G-band) 




Location bias (Sk)g 


±«(5^) (10-3 


pixel) 




13.3 


-1.64 ±0.11 


-6.55 ±0.17 


-0.79 ±0.11 


-2.53 ±0.18 


0.89 ±0.13 


5.41 ±0.17 


14.15 


0.13 ±0.16 


3.67 ± 0.33 


2.07 ±0.17 


7.06 ± 0.33 


3.04 ± 0.20 


12.49 ± 0.34 


15.0 


-0.71 ± 0.27 


-0.69 ± 0.48 


0.35 ± 0.25 


-1.16 ±0.36 


3.42 ± 0.28 


3.94 ± 0.36 


15.875 


-1.18 ± 0.32 


-7.15 ±0.49 


-3.61 ± 0.37 


-13.56 ±0.45 


-3.60 ± 0.43 


-18.86 ±0.48 


16.75 


1.21 ±0.50 


22.88 ± 1.14 


-1.04 ±0.57 


11.65 ± 1.02 


-2.98 ± 0.63 


7.38 ± 1.04 


17.625 


-3.11 ±0.76 


5.06 ± 0.96 


-1.23 ±0.85 


-1.95 ±0.99 


-9.33 ± 0.96 


-3.82 ± 1.19 


18.5 


8.54 ± 1.09 


-4.25 ± 1.24 


2.16 ± 1.25 


4.73 ± 1.50 


9.13 ± 1.54 


-10.50 ± 1.60 


19.25 


2.12 ± 1.62 


35.19 ± 1.95 


1.65 ± 1.82 


-3.24 ± 2.08 


-9.38 ± 2.03 


23.92 ± 2.20 


20.0 


23.62 ± 2.39 


25.12 ±2.89 


5.06 ± 2.80 


17.92 ± 3.27 


20.99 ± 3.50 


1 1.74 ± 3.84 



Table 5. Summary of the measured image location biases in the Gaia operating conditions for different stellar image widths, radiation levels, and different 
levels of mitigation: (i) 'none' corresponds to no mitigation (Section [5.6| , (ii) 'ideal' to the presented forward modelling approach associated to an ideal CDM 
and LSF model and calibration (Section [6.3| , and (iii) 'CDM' to the presented forward modelling approach including the current implementation of CDM 
(Section \6.6) . Note that in the latter case, the optimization of the CDM corresponds to the fully optimized case as described in Sections |6.5| and [6^ This 
optimization was performed for the typical reference image only, the same CDM parameters were used for the two other reference images. 
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K.ererence iniage typ^ 


narrow 


typical 


wide 


Background level (e pixel ^ s ^) 


0.44698 


0.44698 


0.44698 


Readout noise (e ) 


4.35 


4.35 


4.35 


^Vindow size (pixels) 


telemetry 


telemetry 


telemetry 


Trap density (traps pixel 


1 


4 


1 


4 


1 


4 


Mitigation 


none 


none 


none 


Magnitude (G-band) 




Standard errors cr^ 


± Va^ (10^3 pixel) 




13.3 


1.80 ± 0.08 


4.26 ± 0.19 


2.00 ± 0.09 


3.50 ± 0.16 


2.08 ± 0.09 


3.08 ±0.14 


14. 


2.31 ± 0.10 


3.55 ± 0.16 


2.67 ±0.12 


4.56 ± 0.20 


3.16 ± 0.14 


4.43 ± 0.20 


15.0 


3.75 ± 0.17 


10.50 ± 0.47 


4.29 ±0.19 


7.43 ± 0.33 


4.40 ± 0.20 


8.59 ± 0.38 


15.875 


7.10 ± 0.32 


17.58 ±0.79 


6.59 ± 0.29 


12.91 ± 0.58 


7.62 ± 0.34 


11.69 ±0.52 


16.75 


9.02 ± 0.40 


17.78 ±0.80 


9.48 ± 0.42 


14.71 ± 0.66 


10.05 ± 0.45 


15.36 ± 0.69 


17.625 


13.04 ± 0.58 


18.63 ± 0.83 


13.60 ± 0.61 


16.47 ± 0.74 


15.79 ±0.71 


18.95 ± 0.85 


18.5 


16.96 ± 0.76 


20.39 ±0.91 


20.49 ± 0.92 


23.48 ± 1.05 


24.99 ± 1.12 


23.72 ± 1.06 


19.25 


26.44 ± 1.18 


30.43 ± 1.36 


29.02 ± 1.30 


32.12 ± 1.44 


33.86 ± 1.51 


34.77 ± 1.56 


20.0 


37.29 ± 1.67 


43.24 ± 1.93 


45.83 ± 2.05 


49.85 ± 2.23 


56.45 ± 2.52 


58.65 ± 2.62 


Mitigation 


ideal 


ideal 


ideal 


Magnitude (G-band) 




Standard errors CTk ± Va,^ (10 ^ pixel) 




13.3 


1.57 ±0.07 


1.91 ± 0.09 


1.77 ±0.08 


2.01 ± 0.09 


1.87 ±0.08 


2.05 ± 0.09 


14. 


2.24 ± 0.10 


2.55 ±0.11 


2.49 ±0.11 


3.09 ±0.14 


3.13 ±0.14 


3.12 ± 0.14 


15.0 


3.45 ± 0.15 


5.90 ± 0.26 


3.90 ± 0.17 


5.01 ± 0.22 


5.18 ±0.23 


4.37 ± 0.20 


15.875 


5.70 ± 0.25 


8.20 ± 0.37 


5.98 ± 0.27 


7.49 ± 0.34 


7.49 ± 0.34 


7.10 ± 0.32 


16.75 


7.92 ± 0.35 


9.72 ± 0.43 


9.15 ± 0.41 


1 1.04 ± 0.49 


9.76 ± 0.44 


10.28 ± 0.46 


17.625 


12.19 ±0.55 


13.36 ± 0.60 


13.51 ± 0.60 


17.43 ± 0.78 


14.26 ± 0.64 


15.39 ± 0.69 


18.5 


16.89 ± 0.76 


18.50 ± 0.83 


19.86 ± 0.89 


23.80 ± 1.06 


22.42 ± 1.00 


24.94 ± 1.12 


19.25 


26.21 ± 1.17 


30.35 ± 1.36 


29.43 ± 1.32 


34.44 ± 1.54 


30.76 ± 1.38 


34.00 ± 1.52 


20.0 


38.64 ± 1.73 


45.00 ± 2.01 


46.82 ± 2.09 


56.58 ± 2.53 


50.36 ± 2.25 


55.76 ± 2.49 


Mitigation 


CDM 


CDM 


CDM 


Magnitude (G-band) 




Standard errors ctk 


± va,, (10^^ pixel) 




13.3 


1.67 ±0.07 


2.72 ±0.12 


1.77 ±0.08 


2.79 ±0.12 


2.03 ± 0.09 


2.65 ± 0.12 


14.15 


2.47 ±0.11 


5.21 ± 0.23 


2.72 ± 0.12 


5.15 ± 0.23 


3.18 ±0.14 


5.38 ± 0.24 


15.0 


4.33 ± 0.19 


7.54 ± 0.34 


3.92 ± 0.18 


5.69 ± 0.25 


4.42 ± 0.20 


5.71 ± 0.26 


15.875 


5.02 ± 0.22 


7.74 ± 0.35 


5.86 ± 0.26 


7.10 ± 0.32 


6.77 ± 0.30 


7.55 ± 0.34 


16.75 


7.91 ± 0.35 


17.99 ± 0.80 


8.95 ± 0.40 


16.05 ± 0.72 


9.96 ± 0.45 


16.51 ± 0.74 


17.625 


12.02 ± 0.54 


15.12 ± 0.68 


13.30 ± 0.59 


15.69 ± 0.70 


15.11 ± 0.68 


18.74 ± 0.84 


18.5 


17.23 ± 0.77 


19.66 ± 0.88 


20.12 ± 0.90 


23.65 ± 1.06 


24.31 ± 1.09 


25.32 ± 1.13 


19.25 


25.55 ± 1.14 


30.77 ± 1.38 


28.33 ± 1.27 


32.95 ± 1.47 


32.16 ± 1.44 


34.73 ± 1.55 


20.0 


37.83 ± 1.69 


45.73 ± 2.05 


45.74 ± 2.05 


51.67 ±2.31 


55.35 ± 2.48 


60.64 ± 2.71 



Table 6. Summary of the measured image location estimator standard errors in the Gaia operating conditions for different stellar image widths, radiation 
levels, and different levels of mitigation: (i) 'none', (ii) 'ideal', and (iii) 'CDM' (see Table|5). 



